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1.0  INTRODUCTION 


Molecular  gas  lasers  have  been  unsurpassed  for  obtaining  high  efficiencies 

with  high  average  power  operation.  For  example,  up  to  30%  has  been 

realized  for  the  CO2  laser,  and  since  this  is  basically  a  three-level  system, 

the  upper  limit  is  the  40%  quantum  efficiency.  However,  diatomic  molecular 

lasers  can  offer  significant  advantages  for  obtaining  operation  at  even  higher 

efficiencies,  and  CO,  a  particularly  promising  candidate,  has  recently 

1-20 

attracted  the  interest  of  several  workers.  Since  a  diatomic  molecule 

has  only  one  ladder  of  energy  levels  corresponding  to  a  single  vibrational 
mode  (in  contrast  to  CO2'  which  has  three  different  modes)  it  is  possible 
to  attain  quantum  conversion  efficiencies  of  nearly  100%.  The  reason  for 
this  is  that  the  lower  laser  level  of  a  given  transition  can  serve  as  the  upper 
level  of  a  subsequent  transition,  thus  permitting  the  vibrational  energy  to 
be  extracted  as  coherent  radiation  through  several  pairs  of  levels  for  which 
the  populations  are  sufficiently  inverted.  For  CO,  it  should  be  noted  that, 
as  a  result  of  the  rapid  redistribution  of  vibrational  energy  among  the  levels, 
a  molecule  reaching  a  terminal  laser  level  can  be  re-excited,  and  thus  the 
absence  of  lasing  on  the  lower  levels  need  not  limit  the  efficiency.  This 
advantage  is  in  contrast  to  the  CO2  system,  for  which  the  energy  of  the 
lower  laser  level  is  rapidly  dissipated  as  heat  energy  through  VT  collisional 
relaxation. 

A  very  serious  problem  which  exists  for  some  diatomic  molecules  (e.g. , 

HF,  HBr,  NO)  is  the  fast  vibration-translation  (VT)  collisional  relaxation 
processes,  which  cause  molecules  to  lose  vibrational  energy  quickly  as  heat. 
Among  the  diatomic  molecules  which  are  useful  for  laser  applications, 

CO  is  unique  in  that  respect,  since  its  VT  lifetimes  are  relatively  long. 

For  a  diatomic  molecule,  the  ladder  of  vibrational  levels  has  energy 
spacings  which  are  nearly  resonant  except  for  the  anharmonic  defect,  and 


thus,  efficient  vibrational  cross -relaxation  processes  are  possible.  In 
fact,  for  the  CO  laser,  it  is  these  rapid  near-resonant  vibration-vibration 
exchange  collisions  (VV)  that  dominate  the  molecular  kinetics  and  pre¬ 
dominantly  determine  the  structure  of  the  vibrational  population 
distribution. 

A  thorough  understanding  of  the  characteristics  of  the  electrically  excited 
CO  laser  involves  molecular  spectroscopy,  radiative  properties,  molecular 
kinetics,  electron  kinetics,  and  discharge  chemistry.  The  general  purpose 
of  this  report  is  to  summarize  our  understanding  of  the  physical  mechanisms 
of  the  CO  laser,  as  well  as  to  provide  a  useful  compilation  of  the  mathematical 
relations  that  are  frequently  required  for  CO  laser  system  calculations. 

These  results  form  the  basis  for  the  molecular  kinetic  model  and  computer 
program  that  was  developed  to  analyze  the  performance  of  an  electrically 
excited,  steady  state  CO  laser  amplifier  system,  and  which  will  be 
described  in  more  detail  in  a  subsequent  section.  A  summary  of  several 
detailed  experimental  measurements  performed  at  NCL  to  determine  CO 
laser  characteristics  under  different  operating  conditions  are  described 
in  a  separate  report. 

The  radiative  characteristics  for  CO  are  treated  in  Section  2.  0,  where 
molecular  spectroscopy,  spontaneous  and  stimulated  emission,  and  laser 
gain  and  saturation  are  discussed.  Section  3.  0  is  devoted  to  a  description  of 
the  pumping  and  relaxation  processes  which  are  important  for  the  molecular 
kinetics,  and  contains  results  of  the  computer  calculations  which  were  per¬ 
formed  for  a  variety  of  parametric  cases.  In  Section  4.  0,  the  electron 
kinetics  of  a  weakly  ionized  plasma  are  treated  briefly,  although  we  have 


2 


not  yet  attempted  numerical  calculations  of  electron  energy  distribution 
functions.  Since  the  plasma  and  molecular  properties  are  coupled 
through  electron-molecule  collisions,  calculations  for  the  electron 
energy  distribution  would  be  desirable  in  order  to  obtain  a  more 
complete  theoretical  understanding  of  electrically  excited  CO  laser 
systems. 
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2.0 


RADIATIVE  AND  SPECTROSCOPIC  CHARACTERISTICS 


2. 1  Energy  Levels  and  Spectroscopy  of  the  CO  Molecule.  The  spectroscopy 

of  carbon  monoxide  has  been  studied  extensively  because  of  its  natural 

occurrence  as  a  combustion  product  in  flames.  An  excellent  up-to-date 
22 

review  can  be  found  in  the  NBS  monograph  entitled  "The  Band  Spectrum 

of  Carbon  Monoxide.  "  Additional  discussions  and  description  of  earlier 

23 

work  can  be  found  in  Herzberg's  classic  book  on  diatomic  molecules. 

Another  useful  reference  for  molecular  spectroscopy  is  the  book  by 
24 

Penner. 

The  energy  level  diagram  of  CO  molecules  is  shown  in  Figure  2. 1.  The 
visible  spectrum  of  CO,  arising  out  of  electronic  transitions  between  the 
upper  electronic  states,  are  dominated  by  three  prominent  band  systems: 


A1*  -  x*  £+ 

Fourth  Positive  System  (2800-1 140A) 

b3X+  -  a3  n 

Third  Positive  System  (2800-1 140A) 

B 1  L+  -  A  1n 

o 

Angstrom  System  (6600-4100A) 

Among  these,  the  A-X  Fourth  Positive  System  in  the  UV  is  the  most 
pronounced  and  is  easily  excitable. 

Although  CO  laser  action  with  low  power  output  has  been  reported  by 

20 

Mathias  and  Parker  on  the  Angstrom  band  electronic  transitions,  the 
primary  interest  lies  in  the  high  output  power  infrared  CO  laser  involving 
transitions  between  the  vibrational-rotational  levels  of  the  ground  electronic 
state.  These  levels  are  described  below  in  greater  detail. 


i 
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Figure  2. 1.  Energy  Level  Diagram  of  the  CO  Molecule 


2. 1.  1  Vibrational -Rotational  Levels.  Figure  2.2  shows  the  vibrational- 

rotational  levels  of  the  ground  electronic  states  of  CO.  Basically,  the 

vibrational-rotational  levels  in  a  given  electronic  state  can  be  considered 

to  arise  from  a  superposition  of  harmonic  oscillator  and  rigid  rotator 

energy  levels,  and  can  be  labeled  with  three  quantum  numbers,  v,  J,  and 

M.  The  vibrational  energy  of  a  harmonic  oscillator  consists  of  the  terms 

G(v)  =  (v  +  1/2)  u  ,  where  u  is  the  vibrational  frequency,  and  v  =  0,  1, 
e  e 

2, . . .  labels  the  vibrational  quantum  levels.  The  energy  of  a  rigid  rotator 

is  given  by  F(J)  =  B  J(J  +1),  where  B,  the  rotational  constant,  is  related 

to  the  inverse  of  the  moment  of  inertia,  and  J,  the  rotational  angular 

momentum  quantum  number,  assumes  the  values  J  =  0,  1,  2,  ...  Each  of  the 

rotational  levels  are  degenerate,  and  are  labeled  by  a  quantum  number 

M  =  -J,  -(J  -  1),  ...  0,  1,  2,  ...  (J  -  1),  J.  The  total  degeneracy  of  a 

rotational  state  with  rotational  quantum  number  J  is,  therefore,  g  =  (2J  +  1). 

J 

This  leads  to  a  spectrum  of  molecular  energy  levels  given  by 

E(v,  J,  M)  =  E(v,  J)  =  G(v)  +  F(J).  (2.  1) 

F  jwever,  in  practice,  the  vibrational  energy  of  molecules  is  more 

accurately  described  by  an  anharmonic  oscillator,  and  (due  to  centrifugal 

forces)  the  rotational  spectrum  is  actually  that  of  a  non-rigid  rotator.  Con- 

23 

sequently,  it  can  be  shown  that  the  vibrational  and  rotational  terms  are 

G(v  )=  (v  +  1/2)  (j  -  (v  +  1/2)2  u  x 
e  e  e 

+  (v+  1/2)3  u  V  +  (v  +  1/2)4*.*  (2.2) 

6  C  v  C 

and 

F(J)  =  BJ(J  +  1)  -  DJ2  (J  +  l)2,  (2.  3) 
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where  D  is  positive  and  «  B,  and  u  z  « u  y  «  <J  x  u  .  Further- 

e  c  coco  o 

more,  it  is  usually  not  possible  to  completely  ignore  the  interaction 

between  the  vibrational  and  rotational  motions.  If  these  two  systems  are 

coupled,  the  energy  cannot  be  given  simply  as  a  sum  of  G(v)  and  F(J). 

Since  the  molecule  is  simultaneously  vibrating  and  rotating,  the  internudear 

distance  is  changing,  and  thus  the  moment  of  inertia,  which  is  related  to  1/B, 

is  also  changing.  However,  because  the  rotational  frequencies  are  low 

compared  to  the  vibrational  frequencies,  it  is  usually  plausible  to  use  mean 

values  B  and  D  for  the  rotational  constants  when  the  molecule  is  in  a 
v  v 

given  vibrational  state  v.  Thus,  if  the  rotational-vibrational  coupling  is  not 

23 

too  strong,  it  is  satisfactory  to  retain  the  expression  (2.1)  for  E(v,  J), 
with  B^  and  given  by 

Bv  *  B^  -  a  (v  +  1/2)  +  y(v  +  1/2)2  (2.  4) 

D  =  D  +  0(v  +  1/2).  (2.5) 

V  C 

For  CO,  the  vibrational  constants  (j  ,  x  ,  u  y  ,  and  u  z  have  been 

.  _  e  e  e  e  e  e  e 

determined  by  Patel  ’  ,  and  the  rotational  constants  a  ,  fi,  V,  B  ,  and 

25  * 

D  have  been  given  by  Benedict  et  al.  These  are  summarized  in  Table  2. 1 

O 

below. 
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we 

2169.  82 

CJ  X 

e  e 

13.  292 

u  y 

e  e 

. 01082 

it)  z 
e  e 

5.  72  x  10‘5 

B 

e 

1.93141 

D 

e 

6.  18  x  10"6 

a 

.  017520 

fi 

-1.  76  x  10“9 

Y 

2.  96  x  10"6 

Table  2. 1.  Vibrational  and  Rotational  constants  of  the 
CO  molecule  in  wavenumbers  (cm-*). 

2. 1.2  Infrared  Transitions  for  CO.  The  first  order  allowed  infrared 
spectrum  for  CO  results  from  transitions  between  the  vibrational-rotational 
levels,  subject  to  certain  selection  rules.  The  selection  rules  for  electric 
dipole  radiation  lead  to  allowed  transitions  for  v'  «  v  +  1,  and  J'  =  J  £  1. 

The  reason  for  this  is  that  the  first  order  electric  dipole  moment  is  pro¬ 
portional  to  the  displacement,  and  is  a  vector.  This  leads  to  the  vector  addi¬ 
tion  rule  for  the  change  in  J  values,  and  to  4  v  =  +1,  characteristic  of  harmonic 
oscillators.  (The  fact  that  the  vibrational  motion  is  slightly  anharmonic  will 
relax  this  selection  rule  somewhat,  although  transitions  for  which 
A  v  t  +1  will  be  much  weaker. )  The  transitions  for  which  (J  -  l)-w  J  or  (J  +  l)-»  J 
are  referred  to  as  P(J)  or  R(J)  branches,  respectively,  with  the  label  J  denoting 
the  angular  momentum  of  the  final  state.  (There  is  no  Q(J)  branch,  correspond¬ 
ing  to  transitions  J,  for  CO.  This  is  a  consequence  of  the  fact  that  the  total 
angular  momentum  is  equal  to  the  rotational  angular  momentum,  with  no 
electronic  contribution.  The  selection  rule  on  J  then  arises  from  a  Clebsch- 


\ 
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Gordon  coefficient,  C(J,  1,  J* ;  M,  . which  vanishes  for  J  =  J'. ) 

The  first-order  allowed  transition  frequencies  for  the  two  branches  are 
given  by 

P(J)  Branch: 

V  =  E(v,  J-l)  -  E(v-1,  J)  =  G(v)  -  G(v-l)  +  I  v(J-l)  -  Fyl(J)  (2.6) 
R(J)  Branch: 

E(v,  J+l)  -  E(v-1 ,  J)  =  G(v)  -  G(v-l)  +  Fy(J+l)  -  Fyl(J)  (2.  7) 

and  these  can  be  combined  into  a  single  expression, 

»=  V  +  (B  +  B  .)  m  +  (B  -  B  .  -  D  +  D  .)  m2  (2.8) 

O  V  V-l  V  V-l  V  v-l 

-2(D  +  D  )  -  (D  -  D  .)  m*„ 

'  v  v-l  v  v-l' 

where  jT  =  G(v)  -  G(v-l),  and  where  m  =  -J  or  (J  +  1)  for  the  P(J)  or  R(J) 
o 

branch,  respectively.  (Throughout,  we  shall  use  the  convention  that  V 
represents  frequency  in  sec  *,  and  V  *  V  fc  represents  frequency  in  wave- 
numbers,  cm**,)  A  computer  program  was  written  to  calculate  the  P(J) 
and  R(J)  transition  frequencies  and  air  wavelengths  for  v  =  1  to  30  and 
J  =  1  to  50,  and  to  sort  them  into  numerical  order.  This  enabled  the 
experimentally  observed  CO  laser  output  lines  to  be  quickly  and  readily 
identified. 


2. 2  Spontaneous  Emission.  From  the  quantum  theory  of  radiation,  the 
transition  rate  for  spontaneous  radiation  of  a  molecule  from  an  initial 
vibration -rotation  state  (v\  J',  M')  to  a  final  state  (v,  J,  M)  is  given 
by  the  Einstein  A-coefficient, 


A  =  64 *>,(vl»  J ».  1  «(v»  j'  M'-»  v  J  M)  ”  (2  9) 

v\  J.M  ,.3  |/|(v,j,m  v,j,M)  \c.v) 
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where  |  ft  (v1,  J*,  M*-^v,  J,  M)j2  is  the  square  of  the  dipole  matrix  element, 
and  where  the  frequency  |/(v',  J1,  M'-*-v,  J,  M),  which  is  independent  of  the 
quantum  numbers  M\  M,  was  given  by  Eq.  (2.  8)  for  v  =  v*  -  1  for  both  P 
and  R  branches.  Since  a  priori  knowledge  of  the  initial  and  final  degenerate 
quantum  numbers  M'  and  M  is  not  possible,  the  total  spontaneous  radiation 
rate  for  the  transition  v'-^  v,  J'-»-  J  can  be  obtained  from  Eq.  (2.  9)  by  averaging 
over  the  initial  M1  states,  and  summing  over  the  final  M  states.  Since  the 
initial  M'  states  all  have  equal  probability  of  occurence,  viz.  1/gj,,  where 
gji  =  (2J'  +  1)  is  the  rotational  degeneracy  of  the  initial  state,  we  obtain 


Z  (1/gji)  Av.t  j«,  M'-^v,  J,M  * 
M'M 


For  CO,  it  can  be  shown  that  (cf.  Ref.  24  ,  p.  130-133), 


(2.10) 


M?M  |<v',J',M'|  U  |  v,  J,  M  >|2  -  |<V  |R|  v  >|2  X 

j  Ji(J,  J'+l)  +  J'd(J+l,  J')j,  (2.11) 

and  from  this  expression,  it  is  apparent  that  the  two  possible  decay  modes 
correspond  to  a  P  and  R  transition.  Thus,  the  transition  rate  (2. 10)  becomes, 
finally, 


64>t4|/3(v«,  J«^.  y,  J) 
3hc3  (2 J'  +  1) 


<  v* 


X 


j  J  i  ( J ,  J 1  +1)  +  J 1 J  ( J+l,  J ')  ,  (2.12) 
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The  total  spontaneous  decay  rate  for  v'-e*v  can  be  obtained  by  averaging  over 
all  initial  J1  levels,  and  summing  over  all  final  J  levels.  Thus,  if  P(v',  J') 
denotes  the  probability  that  the  initial  rotational  level  is  J',  we  obtain 


J'J  P(V’’J')AV'J'-*V,J 


(2.13) 


For  example,  if  the  rotational  levels  are  in  thermal  equilibrium  at  some 

temperature  T  ,  the  initial  probability  coefficients  p  (v\  J')  will  be 
rot 

Maxwell -Boltzmann  factors, 


P(V.J')  -  .xp  +D/kTrot], 

V* 

V1 


(2.  14) 


where 


Q'v°‘*  J.  <2J'  +  *>  “P  [-Bv.J'<J'+1>/kTrot] 


(2.  IS) 


is  the  rotational  partition  function  of  the  initial  vibrational  level  v1.  If  the 
sum  over  J  in  Eq.  (2.  13)  is  carried  out  first,  we  obtain 


|<  v*  1  R  I  v>|  2  £  (P(v',J')/gjt)  X 

J' 

j(J'+l)U3  (v'.j’^v,  J'+l)  +J'  y3(v\  J'*v,  J'-l)  J  ,  (2.16) 


which  involves  the  P(J'+1)  and  R(J'-l)  frequencies.  From  Eq.  (2.  8),  these 
are  givon  by 


|/(v',  J'-*.v,  J'+l)  -  V  (v',v)  -  (J'+1)(B  +  B  ) 

O  V  V 


(2.17) 


V(v\  J'-*-v,  J'-l)  s  |/o(v'  v)  +  J’(Bv,  +  By) 


(2.  18) 
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and  thus,  the  expression  in  brackets  in  Eq.  (2.  16)  becomes 


(J>  +  1)  v\  -  3(J«  +  l)2  U2  (Byt  +  Bv)  +  j'Tf2  + 

o 


3j'2  T2(B  ,+B  )  +  ...a^J'+l)!?3  -  3|?(B  ,+B  )1.(2.  19) 
o  v’  V  L  O  o  v’  V  J 


Since  the  rotational  constants  B^,  and  B^  are  negligible  compared  with 
the  second  term  in  expression  (2. 19)  can  be  neglected,  and  Eq.  (2. 16) 
reduces,  finally,  to 


64rr4  Va  ( v',v )3 


Ir  ,  I2* 

v  'v  9 


(2.  20) 


where  p  (v't  v)  is  the  energy  (frequency)  difference  in  cm  between  the 

levels  v'  and  v,  and  T  ,  denotes  the  spontaneous  lifetime  for  the  radiative 

v',v 

decay  from  v'^v.  (Note  that,  after  the  sum  over  J  in  Eq.  (2.  13)  was  carried 
out,  the  sum  over  J'  collapsed  to  the  trivial  sum  of  p(v',  J'),  which  is  just 
equal  to  unity.  Thus,  the  total  spontaneous  radiative  decay  rate  of  a  level  v' 
is  independent  of  the  distribution  of  the  rotational  J*  levels. )  For  a  harmonic 

i  .  2 

oscillator,  the  vibrational  matrix  elements  |R  .  |  are  approximately 
2 

equal  to  v  | Rjq|  •  For  CO,  it  has  been  shown  that  a  Morse  potential 
describes  the  matrix  elements  |R^  ^  ^  |2  more  accurately  than  a  quadratic 
potential,  and  that  gives  the  result 

|R  ,  /»  I*  =  A-  (l/*.-2v-l) -  x 

I  v'  v  101  v!  v'-v-l  2 

77  (1/  *  -  v'  +  t)  (v*  -  v)  (l/*e-3) 

t=0 


(l/xe  -  2)2/(l/xe  -  v»-v  -  l)2, 


(2.21) 
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where  x  is  a  measure  of  the  anharmonicity  of  the  molecule,  and  is 
© 

defined  by 


x  =  u  x  /  u  (2. 22) 

e  e  e  e 

The  expression  is  valid  for  v1  >v;  for  v1  <  v,  the  corresponding  expression 
can  be  obtained  simply  by  interchanging  v'  and  v  in  Eq.  (2.  21).  It  is  easy 
to  verify  that  this  result  reduces,  to  a  first  approximation,  to  the  form 
|R  .I2  =  v  |R,  J2  which  is  characteristic  of  a  harmonic  oscillator. 

From  Eq.  (2.  20),  it  follows  that  the  transition  rate  A^,  can  be  expressed 
in  terms  of  A^,  and  the  ratio  of  the  squares  of  the  transition  dipole 
moments,  as  given  by  Eq.  (2.21): 


V  *10  <  9.  (v|v> /rio  >3  lRv',v/Riol2  •  <2‘23> 

This  equation,  in  conjunction  with  Eq.  (2.21),  is  useful,  since  it  allows 
the  radiative  decay  rates  for  Av  =  1  and  A  v  =  2  to  be  calculated  for  high 
vibrational  levels  interms  of  the  basic  rate  A^q.  For  these  cases,  Eq.  (2.23) 
becomes 


lRv.v-l/Rl/  =  vt(l/*e  -2)/<l/x.  ■  2v>i2  X 

J  (1/Xe  -  2v  -  l)(l/xe  -  2v  +  1)  | 
j  (l/xe  -  3)(l/xe  -  v)  ) 

and 

lRy.v  -  2/R10*2  =  4“  v(v-11  ((1/xe  -  ‘2v  +  »]2  X 

j  (l/xe  -  2v  -  l)(l/xe  -  2v  +  3)  | 
j  (l/xe  -  3)(l/xe  -  v)(l/xe-v  +1)}  * 


(2.24) 


(2.25) 
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The  value  of  the  transition  dipole  moment  R..  for  the  CO  molecule  is 
,  -20  10 

R^(CO)  =  9.  906  x  10  esu,  which  corresponds  to  a  spontaneous  radiation 
rate 


*10  ■  l/- 


64  jr4  V. 


10 


sp 


3h 


|R1q|2  =  30.3  sec-1 


(2.26) 


2.3  Stimulated  Emission.  In  the  presence  of  external  radiation,  the  rates 
for  stimulated  emission  or  absorption  of  electromagnetic  energy  by  excited 
molecules  will  be  expressed  in  terms  of  the  Einstein  B-coefficients.  The 
probability  per  second  that  incident  monochromatic  radiation  at  frequency  v 
induces  a  molecular  transition  from  an  initial  state  (v\  J'.M1)  to  a  final 
state  (v,  J.M)  is  given  by 


r1  =  Jizl 

v',  J'.M'^v,  J.M  c 


S(v,i/12)b 


v'J'M'-*V,  J,M 


(2.27) 


where  I(P)  is  the  radiant  power  per  unit  area,  S(P,  V^)  is  the  normalized 
line  shape  function  for  the  molecular  resonance  centered  at  frequency 
V  =  l/(v',J'-*.v,  J),  and  B  ,  is  the  Einstein  B -coefficient, 

defined  by 


3 

Bv\  J',M'-*v,  J,M  =  AV>,  J',M'.»V,  J,M/(87fh  *12 


(2.28) 


It  is  a  fundamental  characteristic  of  quantum  theory  that  the  transition  rates 
for  each  elementary  process  and  its  inverse  are  equal,  assuming  that 
certain  invariance  conditions  are  satisfied.  This  is,  in  substance,  the  principle 
of  detailed  balancing,  according  to  which  the  rates  for  forward  and  reverse 
microscopic  processes  are  equal  for  a  system  in  statistical  equilibrium. 
Therefore, 


Bv',J',MVv,J,M  =  Bv,  J,M^v\  J\M'  , 


(2.29) 


15 


Again,  it  can  be  argued  that,  since  there  can  be  no  knowledge  of  the  initial 
of  final  degenerate  quantum  numbers  M  and  M1,  it  is  more  relevant  to  con¬ 
sider  only  the  total  stimulated  rates  for  the  processes  (v1,  J1)-*  (v,  J)  and 
(v,  J)-^  (v‘,  J‘).  These  can  again  be  defined,  as  for  the  case  of  the  A- 
coefficient  in  the  previous  section,  as  an  average  over  initial  and  sum  over 
final  M  states: 


B 

v 


i 


,  J’ 


V,  J 


L 

M'M 


^J'^v'J'M'-^v,  J,M 


(2.  30a) 


t 

MM* 


(!/gj)  Bv,  J.M-^v'J'M' 


(2.  30b) 


It  follows  immediately  from  Eq.  (2.  30)  that  the  reduced  Einstein  B-coefficients 
satisfy  the  relation. 


gT  B 
J  v 


,  v',  J'  ”  ®J'  ®v',  J'-^  v, 


J* 


(2.31) 


Furthermore,  they  can  be  expressed  in  terms  of  the  reduced  A-coefficient, 
given  by  Eq.  (2. 12),  as 


B 

v 


'.J'-#. 


V,  J 


v.J/'8,rh?B 


)• 


(2.32) 


A  particular  case,  which  will  be  required  later  for  the  discussion  of  laser 
gain,  are  the  B-coefficients  for  the  P(J)  and  R(J)  branches  of  a  transition 
v^^v  -  1: 


B 


v,  J  -  1  «♦  v-1,  J 


Hit  1  In  I2 

2  1  v,v-r  ’ 


3h  g 


J-l 


and 


B 


„  Z£  (J  +  i)  ,  a  |2 

v,  J+l—v-l.J  2  'v.v-l'  • 

3h  «J  +  1 


(2.  33a) 


(2.33b) 
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As  Eq.  (2.27)  shows,  the  stimulated  emission  properties  of  the  molecular 
system  depend  upon  the  resonance  line  shape  function  S(|/).  Therefore,  before 
proceeding  with  a  discussion  of  gain  and  saturation  characteristics,  it  will 
be  useful  to  summarize  some  basic  information  about  line  broadening 
mechanisms  in  molecular  gas  lasers. 

2.4  Line  Shape.  One  of  the  important  factors  that  determines  the  stimulated 
radiative  characteristics  of  the  laser  gain  medium  is  the  width  and  shape  of  the 
spectral  emission  line.  In  general,  the  radiative  transition  will  not  be 
monochromatic  in  frequency,  but  will  be  broadened  due  to  lifetime  effects, 
collisional  effects,  or  medium  inhomogeneities.  The  probability  density  that 
the  observed  emission  will  be  in  the  frequency  interval  ( v ,  v  +  dv)  will  be 
described  by  a  line  shape  distribution  function  g(V),  which  shall  always  be 
assumed  to  be  normalized  to  unity: 

•» 

J  dvg(u)  =  1.  (2.34) 

In  this  section,  a  brief  discussion  of  the  broadening  mechanisms  important 
for  molecular  gas  lasers  will  be  given,  with  a  summary  of  several  useful 
relations.  In  general,  a  Lorentz  line  shape  results  from  processes  that 
limit  the  lifetime  of  excited  states,  modulate  their  energy  levels,  or  induce 
nonradiative  transitions  between  them;  it  also  occurs  for  collisional  processes 
that  interrupt  the  phase  of  the  radiating  molecules,  and  thereby  modulate  the 
emission  frequency.  Alternatively,  medium  inhomogeneities  may  produce  a 
distribution  of  emission  frequencies,  and  this  usually  results  in  a  Gaussian 
line  shape.  For  molecular  gas  lasers,  this  can  arise  from  Doppler  shift  of 
the  transition  frequencies  of  molecules  whose  velocity  distribution  is  Boltzmann. 
The  nature  of  homogeneous  broadening  is  fundamentally  different  from  that  of 
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inhomogeneous  broadening;  in  the  former  case,  all  of  the  molecules  in  the 
medium  share  the  same  basic  emission  properties,  while  in  the  latter  case,  all 
of  the  molecules  are  different.  This  has  important  consequences  for  saturation 
effects,  both  in  regard  to  gain  characteristics  and  line  shape  distribution.  In 
the  homogeneous  case,  saturating  radiation  near  resonance  can  be  distributed 
to  all  of  the  molecules  of  the  medium,  while  in  the  latter  case,  different 
molecules  behave  independently.  For  molecular  gas  lasers,  natural  line 
broadening  due  to  spontaneous  radiation  or  nonradiative  decay  induced  by 
collisions  is  always  negligible;  the  only  two  mechanisms  that  are  important 
are  Doppler  and  pressure  broadening,  which  shall  be  discussed  in  more 
detail  below. 


2.4.1  Doppler  Broadening.  For  conventional  gas  lasers,  operating  at  total 
pressures  typically  less  than  20  torr,  inhomogeneous  Doppler  broadening 
is  the  dominant  mechanism  that  determines  the  radiative  linewidth.  If  it  is 
assumed  that  the  molecular  velocities  are  distributed  according  to  a 
Maxwell -Boltzmann  function,  then  the  fractional  number  of  molecules  in 
(vx,  vx  +  dvx)  is  given  by 

dn(v  ,  v  +  dv  )/n  =  (M/2irRT)*^  exp  [»Mv^  /2RT]dv  .  (2.  35) 

X  X  X  X  X 

If  the  x-direction  is  chosen  as  the  line  of  sight,  the  observed  emission  fre¬ 
quency  from  molecules  having  a  velocity  component  vx  will  be  Doppler -shifted 
to 

V  =  Vq  (1  +  vx/c),  (2.  36) 

where  v  is  the  center  frequency  of  the  transition.  The  probability  that  the 
observed  emission  lies  in  the  frequency  interval  (  v  ,  V  +  dP)  is  defined  to  be 
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j 


dl/g _(  V ,  V  )  -  dn(v  ,  v  +  dv  )/n,  10  that  g_(l/,  1/  )  can  be  obtained 
D  o  x  x  x  D  o 

directly  from  Eq.  (2.  35)  by  substituting  the  result  (2.  36)  and  dv^  =  dv/v  , 

where  v  -  V  /c  is  the  center  frequency  expressed  in  wavenumbers 
-1  °  ° 

(cm  ).  Thus,  the  line  shape  function  for  Doppler  broadening  is  given  by 


where 


SD^*  *'o)  = 


1 

Av  "Sir 

D 


exp  [  - 


MLlHoL  1 

z  J  * 

4  4 


(2.  37) 


(2.  38) 


is  the  half -width  of  the  distribution,  measured  at  the  1/e  point.  (Note  that 
the  A  V^  defined  here  omits  the  Vln  2  factor  which  often  occurs  in  the 
literature,  and  which  would  correspond  to  the  half -width  at  the  half -power 
points. ) 


2.4.2  Pressure  Broadening.  At  higher  pressures,  collision  broadening 
provides  the  dominant  contribution  to  the  linewidth.  Collisional  processes 
represent  a  homogeneous  broadening  mechanism,  with  a  line  shape  function 


given  by 


26,27 


gp(l/.  l/)  =  U* /*)[(*'-*'  -#v  f+Av2]"1  ,  (2.39) 

f  o  c  o  c  c 

where  6  is  a  pressure -induced  shift  of  the  center  frequency  p^,  and  where 

A  V  i*  the  pressure-broadened  linewidth.  Except  for  extremely  high 
c 

pressures  ( several  atm),  the  frequency  shift  iv  is  negligible,  so  that 

c 

the  pressure -broadened  line  shape  function  can  be  expressed  by 

gpfi'.  V  )  =  (Av  In)  Uv-v  )2  +Av*]ml  ,  (2.40) 

ir  o  c  o  c 
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where 


Av  =£  /(2w)  (2.41) 

c  c 

is  the  half -width  at  half -power  points.  Here,  f  is  the  collision  frequency, 
given  by 

fc=  X  V,CO’X)nX<Vr.l(CO'X,>T-  ,2-42» 

where 

<Vrel  (C°’X)  >T  =  [lSRT/tr  MMcJ1  +  Mx-1)]1/2  (2.43) 

is  the  relative  thermal  velocity  between  CO  and  X  molecules  at  temperature 
T,  R  =  8.  317  x  107  erg/mole/°K  is  the  gas  constant,  M's  are  molecular 

3 

weights,  n^  is  the  number  of  X  molecules /cm  ,  and  0  (CO,  X)  is  the  optical 
collision  cross  section  for  the  pair  of  molecules  (CO,  X).  The  "hard  sphere" 
collision  cross  section, 

0(X,Y)  =  7T(dx  +dy)Z/4  (2.44) 

determined  from  kinetic  experiments  is  usually  a  good  approximation  for  the 
effective  optical  cross  section,  since  every  collision  interrupts  the  phase  of 
a  radiating  molecule.  (In  fact,  even  "softer"  collisions  suffice,  since  the 
optical  cross  section  is  typically  slightly  larger  than  the  hard  sphere  cross 
section. )  Thus,  if  reliable  data  on  the  optical  collision  cross  sections  is 
lacking,  the  "hard  sphere"  diameters  d  ,  which  are  tabulated  in  the  table 
below  for  several  molecules  of  interest,  can  often  be  used. 
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X 

M 

dx(A) 

co2 

44 

4.00 

CO 

28 

3.59 

N2 

28 

3.68 

Ar 

40 

3.42 

He 

4 

2.58 

Table  2.2 


2.4.3  Combined  Doppler  and  Preasure  Broadening.  In  the  pressure  and 

temperature  region  where  the  effects  of  collision  and  Doppler  broadening 

are  comparable,  a  more  complicated  line  shape  function  results.  For  this 

intermediate  case,  the  radiative  line  shape  function  will  reflect  the  combined 

effects  of  both  broadening  mechanisms,  and  will  be  expressed  in  terms  of  a 

convolution  of  the  two  distributions  given  by  Eq.  (2.  37)  and  (2.  40).  The 

relative  number  of  molecules  whose  emission  frequency  lies  in  (p',  p'  +  dp') 

due  to  Doppler  shift  froma  center  transition  frequency  pQ  is  dp'  g^(p',  PQ); 

the  relative  number  of  molecules  whose  emission  frequency  is  shifted  from 

V'  to  the  interval  (p,  P  +  dP)  because  of  collisional  effects  is  dPgp(p,P'). 

Hence,  the  relative  number  of  molecules  whose  emission  frequency  is  in 

(P,  p  +  dP)  for  a  transition  centered  at  Pq,  Doppler -shifted  toP',  will  be 

given  by  a  distribution  function  S(P,  p  ;  P'),  defined  by  dPS(p,  P  ;P')  = 

o  o 

dPg  (p,  p')  dP'  g  (p1,  P  ).  Integration  over  all  of  the  intermediate 

ir  U  O 

frequencies  P'  will  then  produce  the  line  shape  function  for  combined 
pressure  and  Doppler  broadening, 

S(P,  PQ)  =  /dP'  gp(P,  P')  gD(P’,  PQ)  .  (2.45) 
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Since  all  of  the  distribution  functions  discussed  here  depend  only  upon 
the  difference  of  the  frequencies,  Eq.  (2.  45)  can  be  written 


S(v-vo)-  J  dv'  gp(l*'-|,0  1  -V1)  g3  (V1)  i  (2.46) 

—  «e 

which  is  the  familiar  convolution  integral  form.  To  evaluate  this  distribution 
function,  it  is  convenient  to  introduce  the  Fourier  transforms, 

g  (r)  s  J  dyei*'Tg(i/)t  (2.47) 

which  can  be  shown  to  be 

gp(r)  =  exp  [  i  v^T  -4V.M  ]  (2.48a) 


and 


gD(r)  =  exp  [  i  vr  -  TZAV* /A  ].  (2. 48b) 

for  the  pressure  and  Doppler  broadened  line  shapes,  respectively.  The 

convolution  integral  (2.46)  for  S(p)  can  then  be  expressed  as  the  inverse 

x 

Fourier  transform  of  the  product  of  gp(r)  g^fr) 

to 

.  _ 1  f,_  -i  VT-bV  |t| 

S(V)  -  J  dT  e  c  exp 

-  to 


(-r2i  vd2/4).  (2.49) 
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This  integral  can  be  related  to  the  real  part  of  an  error  function  of  complex 
argument.  The  resonance  line  shape  function  at  the  center  frequency,  S(0), 
becomes 


S{0)  =  —  *  -  J  dr  exp  [  -  Av  t-t^Av  2/4] 

7T  O  C  I 


(2.  50) 


7Z -  «*p(x2)  / 


nAii 


(2.51) 


A  v  vrr 

O 


exp  (x  )  erfc  (x), 


(2.52) 


where 


x  =  AvJA if, 
c  u 


(2.53) 


and  erfc(x)  is  the  complementary  error  function, 


erfc( x  )  *  -&  /-«•' 


(2.  54) 


It  is  easily  verified  that  the  line  shape  factor  at  resonance  reduces  to  the 
results  for  pure  Doppler  or  pressure  broadening,  if  the  limit  x«*»0  or 
x*«»  is  taken  in  Eq.  (2.  52): 


1)  Doppler  broadened  limit:  x«^0,  and 


S(O)-** 


A  V7T 


(2.  55a) 


2)  Pressure  broadened  limit:  x^ao^  and 


S(0)-» 


n  A  v 


(2.55b) 


where  the  asymptotic  relation  exp  (x  )  erfc(x)^  l/(x  '/W  )  has  been  invoked. 
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2.5  Laser  Gain  Characteristics.  Theoretical  discussions  of  the  gain 

6.7  28 

coefficient  have  been  given  by  Patel  for  CO  and  C02  lasers.  The  gain 
can  be  calculated  as  follows:  the  net  increase  dl  in  power /unit  area  of  an 
incident  beam  of  frequency  v  travelling  a  thickness  dx  of  an  active  medium, 
neglecting  spontaneous  emission  (which  is  emitted  over  all  solid  angles, 
and  is  a  characteristic  of  the  medium  with  on  without  incident  radiation) 
is  given  by 


di  =  h„d*[N2  r‘ j-Nj  r\2  ] 


(2.  56) 


where  N.,  N,  are  the  number  of  molecules  per  unit  volume  in  states  1  and  2, 

i  *  *  * 

and  /"  ,,  r_.  are  the  probabilities  per  unit  time  for  induced  transitions 

12  Z1 

from  states  1  and  2  with  the  emission  (or  absorption)  of  a  photon  with  energy 
hv  .  Expressions  for  .  r have  been  given  by  Eq.  (2. 27)  in  a  previous 
section  in  terms  of  the  Einstein  B -coefficients  and  B^y  If  the  gain  is 
defined  by 


a  s 


(2.57) 


we  obtain 


a  =  (hu/c)  (N2  B21  -  Nj  B12  ]  S  (v,  w12)  .  (2.  58) 

For  the  vibrational-rotational  transitions  of  the  CO  molecule,  the  gain  for 
(v,  J  +1)~»  (v-1,  J)  (i.  e,  the  P(J)  and  R(J)  branches  of  the  v-»(v  -  1)  transition) 
can  be  obtained  from  Eq.  (2.  33),  and  becomes 


v,  J+l 


gff  »(v  J±l-v-l,J)  g  |R  ,2  x 
-*-v-l,  J  3hc  J'  v,v-l' 


I  Jj^Jil 

I  gJ+l 


vzhL 


8 


S(  v  ), 


(2.  59) 
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where 


f  J,  for  P(J)  transitions 


(J+l),  for  R(J)  transitions 


(2.60) 


The  population  densities  of  the  intial  and  final  vibrational-rotational  levels 

have  been  denoted  by  n  .  Since  the  rotational  energy  level  spacing s  are 

v,  J 

small  compared  to  the  average  kinetic  energies  of  the  molecules  under 
typical  temperature  conditions  (kT  >60  cm”*  for  T  Z  77°K),  molecular 
collisions  will  rapidly  thermalize  the  population  distribution  of  the  rotational 
levels  in  any  given  vibrational  state.  For  COg  lasers,  with  B  ~ .  4  cm”*,  this 
rotational  cross  relaxation  is  very  fast,  and  the  J  levels  can  always  be 
assumed  to  attain  a  thermal  equilibrium  described  by  a  Maxwell -Boltzmann 
distribution  of  the  form  given  in  Eq.  (2.14),  with  T  t  —  T  (the  kinetic 
temperature).  For  CO,  the  value  of  B  is  —  2  cm”*,  which  is  five  times 
larger  than  that  for  CO^;  thus,  for  CO,  rotational  cross  relaxation  would  be 
expected  to  be  slower,  since  larger  amounts  of  rotational  energy  have  to  be 
exchanged  between  colliding  molecules  to  thermalize  the  distribution.  If 
rotational  cross  relaxation  is  extremely  fast  compared  to  other  pumping 
rates  and  to  stimulated  emission  rates,  then  even  though  gain  may  exist 
for  several  transitions  simultaneously,  the  line  with  the  highest  gain  will 
begin  to  oscillate  first,  and  the  other  rotational  levels  will  be  depleted 
through  a  single  transition.  On  the  other  hand,  if  the  rotational  cross 
relaxation  is  slow,  it  is  possible  for  several  rotational  branches  of  a  single 
v*sv»l  transition  to  oscillate  simultaneously.  Since  rotational  energy 
exchange  occurs  by  collisional  processes,  its  rate  is  directly  proportional 
to  the  gas  pressure. 
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It  is  usually  a  good  approximation  to  assume  that  the  populations  n^  .  are 
given  by  the  Boltzmann  distribution  (2.14)f  with  the  rotational  temperature 
equal  to  the  kinetic  temperature  T: 


v,  J 


n  g_ 
v  J 

~rot 


exp  [-ByJ(J  +  l)/kT]  , 


(2.61) 


where  denotes  the  total  population  density  of  level  v,  and  where  the 
rotational  partition  function  Q*0t  for  level  v  is  given  by 


O'*  =  f  gj  exp  l-Bv  J(J  +  l)/kT] 


(2.62) 


/ 


dx  exp  (-B 


x/kT)  » 


(2.63) 


The  latter  approximation  holds  because  By«  kT  for  typical  gas  kinetic  tempera¬ 
tures.  The  Boltzmann  constant  k  =  .6945  cm”1  /°K  shall  always  be  understood 
Jo  be  expressed  in  these  units  throughout  this  work  unless  otherwise  specified. 


Under  the  present  assumptions,  the  expression  (2.  59)  for  the  P  and  R  branch 
gains  for  the  transition  v  -a  v-1  become,  using  Eq.  (2. 61): 


v,  J'«a  v-1,  J 


A.n  ~  (v,  J'-»v-l,  J)  S(  v  )  2 

-  S.  |  R  ./R.J  X 

J  v,  v-1  10 


8jriT103  kT 


|nvBv  exp  [-BvJ«(J*  +  l)/kT] 

-  nvJ  Bv-1  exp  l-By-1J(J  +  l)/kT]  j  .  (2.  64) 
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One  of  the  interesting  characteristics  of  molecular  lasers  is  that  the  P-branch 
transitions  can  show  gain  even  when  the  total  vibrational  populations  between 
the  initial  and  final  levels  are  not  inverted.  This  situation  is  known  as 
partial  inversion.  It  can  be  shown  that  the  R -branch  transitions  can  never 
exhibit  gain  under  partial  inversion  conditions.  However,  if  the  total  vibra¬ 
tional  populations  are  sufficiently  inverted,  some  R -branch  transitions 
can  show  gain,  although  it  is  always  less  than  that  for  P-branch  lines.  Even 
for  the  case  of  partial  inversion,  there  are  generally  several  transitions  that 
have  suitable  gain  for  laser  oscillation.  However,  due  to  strong  competition 
effects  between  transitions,  only  a  small  number  of  lines  may  result  in  the 
output.  Therefore,  for  CO  laser  systems  pumped  by  methods  that  do  not 
produce  high  total  inversions,  it  is  only  the  P-branch  transitions  that  are  of 
interest.  For  the  P-branch,  the  gain  equation  can  be  obtained  from  Eq.  (2.  64) 
by  setting  J'  =  (J-l)  and  Sj  =  J.  The  structure  of  Eq.  (2.64)  shows  that  the 
gain  is  a  sensitive  function  of  the  populations,  temperature,  and  J  values, 
and  consequently,  the  spectral  output  of  the  laser  depends  critically  on  the 
operating  conditions.  Figure  2.  3  shows  the  dependence  of  the  normalized 
gain  for  CO  as  a  function  of  J-branch  at  T  -  100°K  for  v  =  6  -►  v  =  5  transitions, 
and  for  a  variety  of  relative  Inversions. 

Note  here  that  Eq.  (2.64)  is  valid  for  both  small-signal  or  saturated 
gain,  provided  that  the  line  is  homogeneously  broadened,  i.  e.  ,  it  depends 
only  upon  the  populations  of  the  vibrational  levels.  Saturation  effects 
will  enter  naturally  into  Eq.  (2.64),  depending  upon  whether  or  not  rates 
of  stimulated  emission  for  the  system  are  comparable  to  pumping  and 
relaxation  processes,  and  if  so,  they  would  have  to  be  included  in  the 
kinetic  equations  which  determine  the  vibrational  populations.  For  the 
situation  where  the  small-signal  region  is  inhomogeneously  broadened, 
the  gain  equation  for  saturation  effects  is  not  given  simply  by  (2.  64),  but 
requires  an  integration  over  the  inhomogeneous  distribution.  This  is 
discussed  in  more  detail  for  a  simple  two-level  example  in  Section  2.7. 
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2.6  Radiative  Contribution  to  Molecular  Kinetics.  The  effects  of 


spontaneous  and  stimulated  emission  have  to  be  included  in  the  set  of 
kinetic  equations,  which  will  be  discussed  in  more  detail  later,  for  the 
molecular  vibrational  level  populations.  In  this  section,  expressions 
for  the  radiative  excitation  and  decay  rates  that  will  be  required  later  are 
developed.  In  the  following  discussion,  it  shall  be  assumed  that  the  radiative 
processes  include  spontaneous  emission  with  dv  =  1  and  d  v  =  2,  and  stimi  lated 
emission  due  to  an  arbitrary  number  of  external  radiation  fields  of  intensities 
l(v,  .  J)  at  frequencies  V  near  the  resonant  frequencies  yj*  that  correspond 
to  P(J),  v-*v-l  transitions. 


The  contribution  of  spontaneous  radiation  with  dv  =  1  and  2  to  the  vibrational 
population  of  level  r  is  dn^/dt  =  R^\  given  by 


R*P  =  -n  (A  .  +  A 
r  r  '  r  -»r-l  i 


,)  +  n  0A  +  n  .A  .  (2.65) 

r-2  r+2  r+2-*r  r+1  r+l-wr 


where  the  Einstein  A-coefficients  have  been  given  by  Eq.  (2. 23).  Spontaneous 
decay  represents  a  pure  loss  mechanism,  with  only  downward  transitions,  and 
for  the  CO  laser,  it  can  only  represent  an  important  relaxation  mechanism  for 
high  vibrational  levels  in  very  low  pressure  systems.  For  the  case  of  stimu¬ 
lated  radiation  processes,  there  are  both  downward  and  upward  transitions 
(emission  and  absorption),  with  the  rate  dn^/dt  =  R  for  the  vibrational 
level  r  given  by 


R.tim 


r,  J-l 


r-1,  J 


+ 


I 

J 


r-1,  J 


r,  J-l 


J  nr,  J  rr,  J-*r+l,  J-l  +  ?  "r+1, J-l  ^r+l,  J-l-w  r,  J. 


(2.66) 


Thus, 


R«*m 

r 


-n  S  .  +  n  _S  , 
r  r  r-1  r-1  r-1 


-n 

r  r 


S 

r 


r+1  +  nr+l  ®r+l 


r  (2.67) 
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where 


r-*r-l 


S(v  ,  /)  A1a  B  R  ,  2 
10  r  r,  r-1 

SJrhcfkT)^*  R1q 


£  J  exp  J  -BrJ(J-l)/kT  I(vr),  (2.68a) 


r-l-^r 


S("'“rJ)A10Br.l  R,.,-l  2 

8>rt.c  (kT)  R10 


£  J  exp  j-Br  lJ(J+l)/kTj  I(iT)  (2.68b) 


(Note  that,  since  k  =  .6945  is  expressed  in  cm"*/°K,  the  quantity  hck  =  1.  38  x 

10~^  erg/°K.  Thus,  the  intensities  I  which  are  conventionally  expressed 

2  7 

in  units  of  watt/cm  must  be  multiplied  by  a  factor  of  10  to  convert  them  to 

2 

erg/sec/cm  .)  If  there  is  only  a  single  external  radiation  field  present,  with 
frequency  ,  the  expression  for  the  rates  of  stimulated  emission  and 
absorption  for  level  v  become 

S(v  ,  )  A._B  .  .  _ 

S  s  J  R  ./R-.r  exp  -B  J(J-l)/kT  I(u)  (2.69a) 

T-V-1  8irhckT~lfl3  |  v.v.l  10|  I  *  v 


S(u,  u  )  A.n  B 

v  10  v-1 


v-l-e  v  8irhckT 


,*-l/Rlo|Z  exp  |  •Bv-l,(,+1,/kT|  !<%)  ,2-69b| 


J  R  ,/R 

V.V-l 


It  is  often  convenient  to  define  optical  cross  sections  and  <T^(v,  J) 

by  the  definition, 


a2l(v,  J)  I(vv ) 


J  m 

V.*  V-1 


(2.  70a) 


a12(v,J)I(w;  ) 


V-1  -*  V 


(2.  70b) 


\ 


so  that 


<,12<V-I,= 


S(  v,  vJ)  A..B  J 
v  v  10  v 


8"kT  “»,o 

S(  *  »v  1  \  A10Bv 


I  Rv  0I2  e*P  j-BvJ(J-l)/kTj  (2.  71a) 


ii|Rv,v-l/R10|2  ex»  |-Bv.lJ,J+1)/kT|  ,2'71b) 


8>rkT 


^0 


In  terms  of  these  cross  sections,  the  gain  for  the  v  -»v-l,  P(J)  transition  can 
be  written 


av,J-l»v-l,  J  *  \  a21,v',)  ■  Vl  ff12<v'J> 


(2.  72) 


The  cross  section  is  often  a  useful  parameter  for  discussing  the  saturation 
characteristics  of  the  laser  medium.  Eq.  (2.67)  and  (2.  65)  will  be  required 

latpr  for  the  discussion  of  the  kinetic  equations  that  determine  the  vibrational 

-i 

level  populations. 


2.7  Saturation  Effects.  When  the  rates  of  stimulated  emission  and 
absorption  become  comparable  to  or  greater  than  rates  of  pumping  and 
relaxation,  the  laser  medium  will  display  saturation  characteristics, 
since  then  the  level  populations  of  the  material  medium  are  substantially 
affected  by  the  radiative  processes.  In  the  limit  that  external  radiation 
fields  are  small,  stimulated  radiation  can  be  neglected  in  the  kinetic 
equations  that  determine  the  vibrational  population  densities.  For  the 
CO  laser  system,  the  effects  of  saturation  and  cross -saturation  will  be 
discussed  in  more  detail  in  a  subsequent  section.  At  this  point,  a  simple 
derivation  of  the  saturation  characteristics  of  a  homogeneously  and 
inhomogeneously  broadened  two-level  system  will  be  presented  in  order 
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to  summarize  some  of  the  basic  considerations  that  are  involved.  For  the 
low  pressure  regime,  much  of  the  experimental  results  for  the  CO 
laser  have  typically  been  in  the  Doppler  broadened  region,  so  measure¬ 
ments  of  saturation  parameters  have  to  be  interpreted  differently  than 
would  be  the  case  for  the  pressure  broadened  region. 

The  rate  equations  for  the  populations  of  a  simple  two-level  system  will  be 
assumed  to  be  of  the  form 

”2  =  -VT2  +  R2  *  (N2  °2l  *  N1  °UW  »  >/h"  <2'  73a> 

N1  =  "V  T1  +  R1  +  (N2  *21  “  N1  *12)I(  V  )/hu  (2-  73b) 

where  Tj  and  are  relaxation  rates,  R^  and  are  pumping  rates,  and 
and  are  cross  sections  for  stimulated  emission,  analogous  to  those 
defined  earlier  in  Eq.  (2.  70).  We  shall  assume  initially  that  the  system  is 
homogeneously  pressure  broadened,  with  a  Lorentzian  line  shape  function 
gp(u  ,  u^),  where  is  the  resonance  frequency  between  levels  1  and  2.  The 
form  of  gp  was  given  by  Eq.  (2.  40)  for  a  collision  broadened  line.  For  the 
moment,  we  shall  think  of  the  two  levels  as  being  vibrational  levels  of  the  CO 
laser,  although  the  kinetics  for  the  CO  system  are  much  more  complicated 
than  the  simple  model  to  be  discussed  here.  Nevertheless,  the  saturation 
characteristics  can  be  understood  best  by  an  initial  appeal  to  a  simple 
derivation,  which  will  serve  to  illustrate  the  basic  concepts.  In  the  steady 
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state,  the  solution  of  these  equations  for  the  gain,  a  (I)  =  (N„  0j2), 

becomes 


a  (I.  w) 


R2  T2  *21 


R  T  0 
1  1  12 


1+  |JL(T1<,12+T2ff21) 


(2.  74) 


(2.75) 


where  a  ,  the  small-signal  gain  corresponding  to  negligibly  small  external 
o 

field  intensity  I(v),  is  given  by 


a 

o 


=  R2T2ff21 


-RlVl2 


(2.76) 


and  where  the  saturation  intensity  I  (v)  is  defined  by 

8 

I8(v)  =  hw/ ( Tj  (w)  +  r2a21(v)).  (2.  77) 

The  physical  interpretation  of  the  saturation  intensity  is  that  aQIg 
represents  the  maximum  possible  extraction  of  power  per  unit  volume 
from  the  homogeneously  broadened  laser  medium. 


As  Eq.  (2.70)  shows,  the  stimulated  emission  cross  sections  0  and  0^2 
depend  upon  the  frequency  v  through  the  line  shape  factor.  For  a  purely 
homogeneously  broadened  medium,  therefore,  the  gain  will  saturate 
according  to  Eq.  (2.  75)  as  the  input  radiation  increases  in  intensity.  If 
we  assume  that  the  levels  2  and  1  represent  v  and  v-1  in  the  CO  system, 
and  that  the  transition  i>j2  is  the  P(J)  branch,  then  the  cross  sections 
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will  be  given  by  Eq.  (2.  71).  Let  Pj  and  p^  be  defined  aa  the  rotational 
Boltzmann  factors,  Eq.  (2. 14),  so  that  the  cross  sections  can  be  written 

*2  ff21(,)  ~~  *p<*  0  <h»/c)  02  -8VJ  -  lRv>  Tj|2  (2. 78a) 

gjV1  =  *?<”•  <V  >  <h»/c>  <*i  — 81-  -  lRv,v.l|2-  <2-78b> 

Then  the  expression  for  saturated  gain  becomes 

08p(v.  vj) 

a  (I, v)  =  - - - - - j  ’  (2.79) 

1  +  yi(u)gp(u,  v J) 


where 


^  =  (R2T2  P2/gZ  “  R1  T1  V8!*  3hCV  ^Rv,vJ  *  (2,80) 


and 


'  =  W.vJ- 

3n  c 


(2.  81) 


Using  the  collision -broadened  expression  Eq.  (2.40)  for  the  line  shape 
function  g(u,  )  gives 

(0Av  /n  ) 

o  (I,  w)  =  - -j-? - = -  (2.  82) 

(v  -„J)  +  4V  +  ri(v)Av  hr 

V  c  c 

which  shows  that  the  gain  line  is  broadened  (i.  e. ,  its  shape  is  changed)  as 
the  intensity  is  increased.  The  important  thing  to  notice  from  Eq.  (2.82)  is 
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that  the  width  of  the  Lorentzian  depends  upon  the  intensity  I  --  that  is,  in  the 
region  where  saturation  effects  become  important,  the  homogeneous  line  is 
both  collision  and  intensity  broadened.  The  form  of  the  saturated  gain 
equation,  (2.  75),  has  been  derived  here  for  a  simple  two  level  system,  but 
it  is  typical  of  the  behavior  in  many  systems. 

If  there  is  an  inhomogeneous  distribution  of  frequencies,  then  the  saturated 
gain  expression  (2.82)  becomes 

?  0g  (v,  v')  j 

0(1,  v)  s  /  du'  - - -  g  (v\  /).  (2.83) 

-  w  1  +  ?I(v)  gp  (v.  u') 


This  equation  is  just  the  convolution  of  a  Gaussian  distribution  g^,  given 
by  Eq.  (2.  37),  and  a  broadened  Lorentzian,  whose  form  was  given  more 
explicitly  in  Eq.  (2.  82).  As  was  shown  in  Section  2.  4.  3,  the  convolution 
can  be  evaluated  in  terms  of  the  complementary  error  function  for  the  case 
of  exact  resonance,  v  =  •  Using  Eq.  (2.  52),  we  obtain 


a  (I**'  12) 

where  the  parameter 


„  frxp  (x2)  erfc(x) 

(  I  +1/:  )1/Z 

D  8 

x  is  defined  as 


(2.  84) 


x 


in  terms  of  the  effective  homogeneous  width, 


(2.85) 


The  parameter 


(1  +  I/I8). 


(2.  86) 


I,  =  nAvc/Y 


(2.  87) 
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represents  the  saturation  intensity  for  the  homogeneously  broadened  line 
at  exact  resonance,  and  in  the  derivation  of  Eq.  (2.  84),  account  has  been 
taken  of  the  fact  that  Eq.  (2.  82)  represents  a  Lorentzian  with  an  effective 
broadening  A  |/g For  the  extreme  Doppler  broadened  region,  for  which 
x  «.  1,  Eq.  (2.  84)  for  the  saturated  gain  at  resonance  reduces  to 


(fi  fy/nAv^) 

a(I)=  (1  +  I/I  )  I/* 

8 


(2.  88) 


while  for  x  »  1,  which  is  characteristic  of  homogeneous  broadening, 
Eq.  (2. 84)  becomes 


a  (I)  = 


(0/nAvc) 

aTTTrr 


(2.  89) 


An  interesting  point  to  notice  is  the  following:  if  operation  is  initially  in 

the  extreme  Doppler  region,  the  effective  spectral  "hole  width"  that  is 

burned  out  of  the  frequency  distribution  by  a  radiation  field  is  jurt  A  v  ,  and 

c 

the  gain  begins  to  saturate  according  to  the  expression  (2.  88).  However, 

as  the  radiation  field  intensity  I  begins  to  increase,  the  effective  homogeneous 

linewidth  given  by  Eq.  (2.86)  increases,  and  when  the  ratio 

x  =  Aveft/Av.Q  becorne8  comparable  to  or  greater  than  unity,  the  saturation 

behavior  goes  over  to  the  homogeneously  broadened  limit,  characterized  by 

Eq.  (2.89).  Notice  that  the  quantities  (  fi /Jn  )  and  (  fi  /nAv  )  that  occur 

c 

in  the  numerators  of  Eq.  (2.  88)  and  (2.  89)  are  just  equal  to  the  resonance 

small-signal  gains  that  are  characteristic  of  the  Doppler  or  pressure 

broadened  regimes  (since  1  fjn  A  vn  =  gn(°)»  and  l/TMv  =  g-r,(0)» 

J-/  u  c  * 

respectively).  Since  the  saturation  behavior  goes  over  to  the  form  (2.  89)  as 
the  intensity  I  becomes  very  large,  I  retains  the  physical  significance  that 
(X  represents  the  limiting  extractable  energy  per  unit  volume  from  the 
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laser  medium.  Note,  however,  that  to  retain  this  interpretation,  ctQ 
must  be  understood  to  be  the  collision  broadened  small-signal  gain  -- 
i.e. ,  <*o  =  0gp(O )=fi/nAvc. 

Since  the  actual  experimentally  measured  small-signal  gain  is  a  measure 
of  the  combined  line  shape  function  S(0)  given  by  Eq.  (2.  52), 


0Q  *  (0A/5rdvD)  exp  (Av*/  Avd  )  erfc  (4uc/4uD>, 


(2.90) 


it  then  follows  that  for  any  arbitrary  region  of  values  Av  ,  A  v_,  and  I, 

c  D 

that  the  gain  <* (I),  given  by  Eq.  (2.  84),  saturates  at  resonance  according  to 


(2.91) 


which  reduces  to  the  expressions  (2.  88)  and  (2.  89)  in  the  appropriate 
limits.  It  further  follows  that,  if  is  the  actual  small  signal  gain,  then 
the  limiting  extractable  energy  per  unit  volume  is  given  by 

ao  Is(  A  vD^A  vc]  exP  ^  v^/  4  £  )/erfc  (AvJ  A  vD>  •  (2.92) 

If  the  experimental  measurements  of  small-signal  gain  and  saturation  intensity 
are  carried  out  in  the  Doppler  broadened  region,  this  expression  reduces  to 


ao  Is  {Av^/yJir  Avc  )•  (2.93) 

Physically,  the  meaning  of  (2.93)  is  that,  for  the  Doppler  broadened  region, 
only  the  fraction  ( irAv  /y/n  Av-)  of  the  spectral  line  is  being  tapped  by  the 
external  radiation  field.  Of  course,  if  the  experimental  data  corresponds 
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to  the  pressure  broadened  region,  expression  (2.  92)  for  the  extractable 
energy  reduces  immediately  to  the  familiar  form,  a  I  . 


3.  0  MOLECULAR  KINETICS:  PUMPING  AND  RELAXATION  MECHANISMS 


3.1  Introduction.  The  program  initiated  at  NCL  for  the  study  and  develop¬ 
ment  of  low  pressure  electrically  excited  CO  lasers  has  resulted  in  a 
considerable  understanding  of  the  relative  importance  of  the  operative 
mechanisms  and  excitation  processes  that  determine  the  vibrational 
populations  of  CO.  In  support  of  the  experimental  research,  a  theoretical 
analysis  was  undertaken  to  develop  a  steady-state  kinetic  model  of  an 
electrically  excited  CO  laser  system  that  would  make  it  possible  to  predict 
the  CW  performance  to  be  expected  under  arbitrary  experimental  conditions 
and  variation  of  parameters.  Because  of  its  high  efficiency,  demonstrated 
with  conventional  longitudinal  electric  discharge  configurations,  the  CO 
laser  is  a  promising  candidate  for  high  power  electric  discharge  systems 
employing  more  sophisticated  techniques  of  plasma  control  for  optimization 
of  operation.  In  terms  of  excitation  and  relaxation  mechanisms,  gas 
constituents,  and  quasi-continuous  mode  of  operation,  long -pulse  transversely- 
excited  E-beam  devices  that  are  presently  being  explored  for  the  CO  laser 
will  be  fundamentally  only  an  extension  of  the  low  pressure  longitudinal 
devices.  Therefore,  the  construction  of  such  a  kinetic  model  represents 
a  valuable  approach  for  extrapolating  low  pressure  results  to  determine 
design  criteria  and  performance  characteristics  for  future  high  energy 
devices. 

A  description  of  the  steady-state  kinetic  model  for  an  electrically  excited 
CO  laser  amplifier  system,  and  the  computer  program  that  has  been 
developed  for  the  calculation  of  vibrational  populations,  gains,  saturation 
parameters,  rates  of  energy  transfer,  and  quantum  efficiency,  will  be 
given  below,  with  a  discussion  of  the  numerical  techniques  employed.  At 
the  present  time,  the  status  of  this  work  is  relatively  complete,  although 
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further  effort  is  required  to  determine  certain  important  constants  and  to 
make  comparisons  between  calculations  and  experimental  measurements. 
Basic  experiments  for  that  purpose,  such  as  small  signal  gain  measurements 
and  double -probe  and  double  Q- switching  techniques  for  measurement  of 
vibrational  cross  relaxation,  are  discussed  more  fully  in  a  companion  report  . 
Although  there  is  presently  a  need  for  improvement  in  the  knowledge  of  the 
values  and  functional  dependence  of  certain  parameters  o£MU  model,  the 
program  has  been  developed  to  handle  the  most  general  posHLle  input 
operating  conditions.  The  computer  code  has  been  compleMK«debugged , 
and  is  performing  successfully.  In  subsequent  sections,  a  Ascription  of 
the  mechanisms  that  have  been  included  will  be  given,  and  representative 
results  (based  on  preliminary  assumptions  concerning  the  rate  constants) 
will  be  presented. 


3.  2  Steady-State  Kinetic  Model.  As  a  general  theoretical  model  of  the 

CO  laser  medium,  we  consider  a  mixture  of  gases  A,  B,  C,  . . . . ,  excited 

by  an  electric  discharge,  where  A  represents  the  active  diatomic  gas,  CO, 

B  a  second  diatomic  gas,  and  C, . . . .  other  monatomic  constituents.  The 

energy  levels  and  populations  of  A  and  B  will  be  denoted  by  E^,  n^  (r  =  0, 1,  2, . . . ) 

and  F_,  N  (R  =  0,  1,  2, . . . )  respectively.  It  shall  be  assumed  that  the  walls 
K  K 

of  the  discharge  volume  represent  a  heat  reservoir  at  temperature  T  with  an 

infinite  heat  capacity,  and  that  the  distribution  of  translational  energies  of 

each  of  the  molecular  species  is  Boltzmann  at  temperature  T.  Likewise, 

the  rotational  energy  levels  in  each  vibrational  state  are  assumed  to  be 

in  good  thermal  contact  with  the  translational  mode  through  molecular 

collision  processes,  so  that  we  may  assume  T  =  T.  The  population 

distributions  n  and  N  of  these  diatomic  systems  are  determined  by  a 
r  K 

variety  of  competing  pumping  and  relaxation  processes,  and  satisfy  the 
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master  equation, 


dn  /  dt  =  R 
r 

+  R 


*  +  KVrT  (Af,  A)  +  R^T(A^,  B,  C. . . . 

VVT(At,Bt)+R,p+R,tim  +  Rmi’C 
r  r  r  r 


VVT  f  4 
+  R  (A  ,  AT] 


=  0, 


(3.1a) 


with  a  similar  equation  for  N  (R  =  0,  1,  2. . . ).  Although  the  ultimate 
source  of  vibrational  excitation  is  by  electron  collisions,  the  method  of 
energy  input  need  not  be  (and  in  fact,  is  not)  the  dominant  mechanism 
in  determining  the  resulting  vibrational  population  distributions.  In  the 
order  of  their  appearance,  the  terms  occurring  on  the  RHS  qf  the  master 
equation  represent  excitation  and  relaxation  rates  for  the  following 
processes: 


i)  electron  impact  collisions 

ii)  VT  collisions  (single  quantum)  by  impact  with  A 

iii)  VT  collisions  (single  quantum)  by  impact  with  B,  C, .  . . 

iv)  VVT  near  resonant  exchange  collisions  with  another  A*  molecule 

v)  VVT  near  resonant  exchange  collisions  with  a  B  '  molecule 

vi)  spontaneous  decay  (Av  =  1  and  Av  =  2) 

vii)  stimulated  emission  processes,  resulting  from  input 

radiation  for  an  arbitrary  number  of  vibrational -rotational 
line  8. 

viii)  miscellaneous  pumping  or  decay  terms 

4 

In  the  steady  state,  the  time  derivative  can  be  set  equal  to  zero,  yielding 
a  multilevel  set  of  coupled,  nonlinear,  algebraic  simultaneous  equations: 

R*  +  RVrV,  A)  +  R^T(At,  B,  C, .  . . )  +  B?VT{Af,  A1" ) 


VVT  4  4  sd 

+R  (AT,BT)  +  R  P  +  R 


stim  +  Rmisc 
r  r 


(3.1b) 


The  processes  i)  to  v)  will  be  described  in  more  detail  in  the  following 
sections.  Spontaneous  and  stimulated  emission  have  been  discussed 
thoroughly  in  Section  2.  6  on  radiative  characteristics,  with  the  rates 
and  R®tim  given  by  Eq.  2.65  and  Eq.  2.67  respectively. 

3.  3  Electronic  Excitation.  The  mechanism  for  vibrational  excitation  in 
electrical  discharge  CO  and  CO.  laser  systems  is  by  indirect  electron 
impact  collisions,  which  proceed  via  the  formation  of  intermediate  compound 
excited  states  due  to  resonant  electron  scattering.  Although  the  ultimate 
source  of  vibrational  excitation  is  by  indirect  electron  collisions,  the  method 
of  energy  input  need  not  be  (and,  in  fact,  is  not  in  the  case  of  the  CO  laser) 
the  dominant  mechanism  determing  the  resulting  populational  distribution. 


The  scattering  cross  sections  for  vibrational  excitation  of  several  molecules, 
such  as  CO,  N_,  and  CO.,  etc. ,  exhibit  pronounced  peaks  for  incident 
electron  energies  of  typically  a  few  electron  volts.  Shulz  has  measured 

these  cross  sections  experimentally  for  several  molecules,  with  the  results 
for  CO  and  shown  in  Fig.  3.1.  From  this  Figure,  it  can  be  seen  that 
electrons  with  energies  of  1.  0  to  3.  0  eV  are  most  efficient  for  vibrational 
excitation. 


Since  the  vibrational  energy  levels  for  CO  and  are  separated  by  much 
smaller  amounts  (-—.2  eV),  it  follows  that  the  collision  process  does  not 
correspond  to  a  direct  electronic  excitation.  The  existence  of  the  phenomenon 
suggests  that  the  excitation  can  be  interpreted  in  terms  of  indirect  resonance 
electron  scattering,  which  arises  from  the  formation  of  a  temporary,  short¬ 
lived,  intermediate  compound  state  CO~  or  N  ".  The  existence  of  well- 

Ca 

defined  peaks  in  the  vibrational  cross  sections  can  be  attributed  to  resonances 
of  the  incident  electron  wave  in  the  potential  well  of  the  molecule.  In  this 


42 


Cwn  Section,  aUtiu;  inlti 


interaction  region,  the  electron  wave  function  at  certain  resonance  energies 

is  distorted  and  enhanced,  and  may  form  a  quasi-bound  state  of  positive 

energy  with  the  molecule,  such  as  CO  or  .  There  is,  however,  a 

finite  probability  of  decay,  characterized  by  the  lifetime  T  of  the  complex, 

with  the  widths  of  the  cross  section  peaks  corresponding  roughly  to  the 

32 

inverse  of  the  lifetime  T  of  the  excited  state.  Herzenberg  and  Mandl 
33 

and  Chen  have  investigated  the  problem  of  vibrational  excitation  by  electron 
collision  theoretically,  and  have  obtained  good  agreement  with  experiments 
based  on  this  model. 

Physically,  the  process  of  vibrational  excitation  can  be  understood  as  follows: 
the  incident  electron  with  a  favorable  resonance  energy  is  trapped  in  the 
molecular  potential  well,  forming  a  negative  compound  state  that  lives  for 
some  lifetime  T .  During  that  time,  the  energy  of  the  electron  "dissolves" 
into  the  available  energy  modes  of  the  complex,  and  when  the  electron  is 
finally  re-emitted,  the  molecule  may  retain  some  of  this  energy  in  the  form 
of  an  excited  vibrational  level.  As  shown  schematically  in  Fig.  3.2  the  energy 
levels  of  the  compound  state  are  not  the  same  as  those  of  the  molecule. 
Electrons  with  incident  energies  of  E  ,  E, ,  E_,  . . .  will  have  an  enhanced 
probability  of  being  captured  to  form  the  excited  compound  state,  and  thus, 
it  is  for  these  energies  that  the  cross  section  will  be  expected  to  display 
resonance  peaks.  Basically,  therefore,  the  excitation  can  be  pictured  as 
occurring  via  a  three  state  process: 

1)  e"  (E)  +  CO^CO" 

2)  CO~  lives  for  time  T,  during  which  time  its  energy  may 
be  redistributed  internally  among  various  energy  modes 

3)  CO’-^  CO(v)  +  e'  (E1 2 3) 
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The  total  rate  of  excitation  dw  (r,  k  -*  s,  k')  into  the  continuum  of  final 
«<  3 

states  dk'/(2Jr)  for  the  scattering  processes  in  which  an  incident  electron 
of  momentum  "hie  excites  a  CO  molecule  from  an  initial  vibrational  level 
r  to  a  final  level  s,  represented  by 

e"(k)  +  CO(r)  -  CO(s)  +  e‘(k'j  (3.2) 


can  be  expressed  in  terms  of  the  scattering  matrix  T  as 

dw(r,  k  -*  s,  k')  =  |t  (k,k*)|2i(E.  -  EJ  dk'/(27T)3, 
n  rs  •  1  f 


(3.3) 


where  and  E^  are  the  initial  and  final  energies, 

E.  =  E  +1i2k2/2m 
l  r 

E  =  E  +-h2k'2/2m 

I  8 


(3.4) 


If  f(k)  denotes  the  electron  distribution  function,  so  that  n  f(k)dk  represents 

© 

the  number  of  electrons  per  unit  volume  in  the  range  (k,  k+dk),  where  n 

e 

is  the  electron  density,  then  the  total  rate  of  vibrational  excitation  from 
r  to  s  is  given  by 


R 


e 


r-»s 


ne(27T/li)  [dkf(k)JdE'/(27T)3  x 

|Trs(k,k')j  $(Er +11  k  /2m  -  Es  -1i  k'  /2m) 


(3.5) 


which  is  the  transition  rate  (3.3)  averaged  over  initial  and  summed  over 
final  electron  states.  The  presence  of  the  ^-function  insures  energy 
conservation  when  the  integral  is  carried  out  over  all  k  and  S'.  The 
microscopic  principle  of  detailed  balancing  for  the  process  (3.2)  is 
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OH*) 


i 


expressed  mathematically  by 


|Trg<£,£')|2  =  |TBr(ic',it)|2,  (3.6) 

so  the  rate  F *  for  the  reverse  process  can  be  written  exactly  the  same 

as  Eq.  (3.  5)  with  f(jc' )  replacing  f (ic).  Eq.  (3.  5)  demonstrates  that  the 

rate  of  vibrational  excitation  depends  upon  the  electron  kinetics  of  the 

plasma,  since  the  electron  distribution  function  f(k)  is  involved.  If  a 

Boltzmann  distribution  is  assumed  for  the  electrons,  characterized  by 

a  temperature  T  , 
e 

f(v)  =  (m/2JTkT  )3/2  exp  (-mv2/2kT  )  (3.7) 

e  e 


where  v  =  hk/m  is  the  electron  velocity  ,  then  it  is  easy  to  verify 
immediately  from  Eq.  (3.  5)  and  the  corresponding  equation  that  would 
apply  for  the  inverse  process  that 


z  E  ~  E  i 
,e  /  s  r  „e 

t  =  exp  l  — — - )  R  , 

r-»s  r  '  kT  '  s~r’ 


(3.8) 


which  is  the  usual  thermodynamic  expression  of  the  principle  of  detailed 
balancing  for  an  electron  gas  that  is  in  thermal  equilibrium  at  temperature  Tg. 


It  is  convenient  to  relate  the  more  familiar  quantities  such  as  cross  section 
and  differential  cross  section  to  the  scattering  matrix  expression  given 
above.  The  scattering  matrix  |t  (lc,  k*)|  is  a  function  of  the  initial 

IT  8 

electron  energy  (and  thus,  k),  and  of  the  final  scattering  angles,  specified 

2 

by  If  the  volume  element  dk  in  Eq.  (3.  5)  is  replaced  by  k1  dk'dflj*,, 

and  the  integral  over  k'  carried  out,  the  differential  scattering  cross 
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flection  can  be  defined  as 


viarJiav  -  T  lTr.<k-flE')|2p,Ef)  «3-9' 

where  P(E^),  the  density  of  final  states,  is  defined  by 

p(Ef)  =  (1/2JT)3  f dk'k'24(Er  +H2k2/2m  -  Ea  -fc2k'2/2m)  (3.10) 

and  v  is  the  initial  electron  velocity.  It  is  easy  to  verify  that 

P(Ef)  =  (m2/h3)v'  (3.11) 

where  v*  is  the  final  velocity  given  by 

v'2  =  v2  +  2(Ef  -  Ea)/m,  (3.12) 

2 

and  that  P( EJ  vanishes  if  2(E  -  E  )/m>v  .  That  is,  if  the  excitation 

i  s  r 

process  corresponds  to  a  vibrational  change  from  r  -*  s,  the  incident 
electron  must  have  a  minimum  energy  of  at  least  (Eg  -  E^)  if  s  >  r. 

For  de-excitation,  of  course,  there  is  no  minimum  energy  requirement 
for  the  incident  electron. 

The  total  cross  section  is  defined  as  an  integration  over  all  solid  angle 
of  the  differential  scattering  cross  section: 

a  (E)  =  /dC(d<7  /dfl)  (3.13) 

r-*s  j  r-*a 
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This  cross  section  corresponds  to  the  experimentally  measured  quantities, 
such  as  those  obtained  by  Schulz  for  CO  and  N^.  which  were  shown  in 
Fig.  3.1.  It  can  be  seen  from  Eqs.  (3.  9)  and  (3.13)  that  the  total  excitation 
cross  sections  for  forward  are  reverse  processes  can  be  related  by 

(v/V)  O  (E)  =  (v'/v)  a  JE‘),  (3.14) 

r-*B  B-»r 

which  is  yet  another  form  for  expressing  the  principle  of  detailed  balancing. 
If  the  Boltzmann  distribution  of  Eq.  (3.7)  is  assumed  for  the  electrons,  it 
can  be  shown  that  the  resulting  distribution  in  energy  will  be  given  by 


and  the 

R 

where 

<v  >  =  (8kT  /»rm)1/2  (3.17) 

e  e 

is  the  average  thermal  electron  velocity  for  temperature  Tq.  It  is 

instructive  to  notice  that  by  using  the  detailed  balance  relation 
2  '2 

v  s  v  <7f  (E1)  'n  conjunction  with  the  energy  conservation 

condition  (3.12),  Eq.  (3.16)  can  be  manipulated  to  demonstrate  the 
thermodynamic  expression  of  detailed  balancing  that  was  given  by 
Eq.  (3.8).  The  net  rate  of  electron  excitation  of  level  r  is 


f(E)  =  — pr—  (kT  )"3/2  E1/Z  exp  (-E/kT  ) 
y/n  e  ® 


(3.15) 


resulting  excitation  rates  R  (s>r)  will  become 

«s 

®  =  n  < v  >  / (kT  )2  f  dEO  (E)  E  exp  (-E/kT  ),  (3.16) 

r-»s  e  e  e  fE  -E  )  r"*“  6 

'  i"  t 
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dn  /dt  =  Re,  where  Re  is  given  by 


R6  = 
r 


-  n  ^  R6  +  ^  n  Re 

r  s  r-»s  s  s  s-»r 

n.  ?[”.  -  nr  'XP^  rkTe'8) 


<v(E)<7  (E  )>_ 
s-*r  T 


(3.18) 


Unfortunately,  the  available  experimental  data  on  electron  impact 
excitation  of  the  molecular  vibrational  levels  is  restricted  to  collision 
processes  for  which  the  molecule  is  initially  in  the  ground  state.  For 
CO  and  N^,  the  data  of  Schulz  includes  cross  sections  for  processes 
0  -» v,  where  v  =  1,  2,  ...  8.  Therefore,  for  the  construction  of  a 
molecular  kinetic  model,  some  assumptions  must  be  made  on  excitation 
rates  for  which  the  initial  molecular  level  is  not  zero.  Based  on  the 
simple  schematic  picture  of  the  resonance  scattering  process  discussed 
above,  and  represented  by  Fig.  3.  2,  it  is  probably  reasonable  to  assume 
that  for  the  first  few  levels,  say  r  =  1,  2,  ...  6,  the  cross  sections  are 
about  the  same,  but  with  the  position  of  the  energy  peaks  shifted  by  the 
initial  energy  of  the  molecular  state: 


*  .(E)  = 

o-*  1 


r-»r  +i 


(E  -  Er) 


(3.19) 


Eq.  (3.16)  can  be  numerically  integrated  to  give  the  electron  pumping 

rates  for  CO.  The  rates  Re  / n  =  <v(E)a  (E)>  obtained  for 

r  -*8  e  r-*s  X 

representative  temperatures  T^  and  various  levels  rf  s  are  shown  in 

Table  3.1.  As  Eq.(3.16)  clearly  demonstrates,  the  electron  excitation 

rates  depend  very  critically  upon  a  favorable  overlap  between  the  electron 

energy  distribution  function  and  the  cross  sections,  a  (E).  For  CO, 

r-»s 

Ng,  and  CO^,  it  is  fortuitous  that  the  cross  sections  are  appreciable  for 
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electron  energies  of  a  few  eV,  which  is  characteristic  of  typical  electric 
discharge  conditions.  (NO,  for  example,  does  not  form  a  compound 
state  NO  ,  and  thus  cannot  be  exploited  as  efficiently  as  CO. )  To  the 
extent  that  it  is  possible,  it  is,  therefore,  always  obviously  desirable 
to  design  the  laser  system  parameters  in  such  a  way  that  this  overlap 
effect  will  be  maximized.  The  optimization  of  the  plasma  characteristics 
in  electric  discharge  lasers  is  an  important  factor  that  determines  the 
efficiency  of  electrical  energy  transfer  into  the  vibrationally  excited 
states  of  the  active  molecules,  and  is,  therefore,  a  fundamental 
consideration  for  increasing  the  conversion  efficiency.  There  are 
several  ways  in  which  the  plasma  characteristics  can  be  altered  to 
enhance  the  laser  operation;  some  of  them  are  relatively  simple,  such 
as  altering  the  discharge  constituents  to  optimize  the  gas  mixture,  while 
others  are  relatively  sophisticated.  Further  discussion  of  some  of  these 
considerations  will  be  given  in  a  subsequent  section  devoted  to  the  electron 
kinetics  of  a  weakly  ionized  plasma. 

3.4  Vibrational  Relaxation  and  Cross -Relaxation.  Collisions  between 
excited  molecules  can  result  in  excitation  or  decay  of  vibrational  quanta 
by  so-called  V-T  processes,  such  as 

A(r+1)  +  A2  A(r)  +  A  +  4E  =  E  ,  -  E 

r+1  r 

B(R+1)  +  B  B(R)+ B  +  4E  =  Fr+1  -  Fr 

A(r+1)  +  B  2  A(r)  +  B  +  dE  =  E  -  E 

r+l  r 

B(R+1)  +  A  =  B(R)  +  A  +AE  =  Fp+1  -  FR 
in  which  the  excess  energy  is  supplied  or  provided  by  kinetic  energy, 
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or  by  means  of  the  so-called  V-V-T  processes,  such  as 

A(r+1)  +  A(s)  **  A(r)  +  A(s+1)  +AE  =  E  -E  +  E  -E  ,  (3,21a) 

'  '  r+1  r  s  s+l 

A(r+1)  +  B(S)  2  A(r)  +  B(S+1)  +AE  =  Ef+1  -Ef  +  Fg  -  Fg+1  (3.21b) 

B(R+1)  +  B(S)  2  B(R)+  B(S+1)  +4E  =  FR+1~FR  +  Fg  -  Fg+1  (3,21c) 

A(r+1)  +  A(s+1)  Z  A(r)  +  A(s)  +AE  =  Ef+1  -Er  +  Efl+1  -Es  (3.21d) 

A(r+1)  +  B(S+1)  2  A(r)  +  B(S)  +AE  =  E  .  -E  +  F_  -  F_  (3.21e) 

IT+l  r  D+l  D 

B(R+1)+ B(S+1)  Z  B(R)  +  B(S)  +AE  =  FR+rFR  +  FS+1  “  FS  (3*21f) 


in  which  vibrational  quanta  of  energy  are  exchanged  between  molecules, 
with  the  excess  energy  transferred  to  or  from  the  kinetic  energy.  The 
rates  for  V-T  processes  decrease  rapidly  as  the  amount  of  kinetic 
energy  4E  that  must  be  absorbed  by  the  translational  mode  increase  &  . 
Furthermore,  the  collision  rates  are  predominantly  determined  from 
dipole  transitions  for  which  dv  =  1,  so  it  is  usually  reasonable  to 
neglect  all  V-T  processes  but  those  which  convert  only  a  single 
vibrational  quantum  of  energy.  Likewise,  V-V-T  transitions  for 
processes  other  than  those  in  which  each  molecule  changes  by  a  single 
quantum  can  be  neglected. 

In  the  first  group  of  V-V  processes,  given  in  Eqs.  (3.21a,  b,  c),  one 
molecule  loses  energy  while  the  second  gains  energy;  thus,  for  those 
cases,  the  excess  energy  dE  exchanged  with  the  translational  mode  will 
be  small.  These  represent  near-resonant  exchange  collisions.  In 
contrast,  the  second  group  of  V-V-T  processes  are  nonresonant  exchange 
collisions,  since  both  molecules  gain  or  lose  energy,  all  of  which  must 
be  transferred  to  or  from  kinetic  energy.  Since  the  vibrational  energies 
of  CO  or  N_  are  ~2000  cm  \  the  translational  energy  that  must  be 
exchanged  is  —4000  cm"  -80  kT  (for  T  -77  K,  for  example).  Thus,  the 
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transition  rate  for  such  processes  as  represented  by  Eqs.  (3.21d,  e,  f) 
would  be  expected  to  be  very  small,  and  they  can  be  neglected.  Only  the 
near-resonant  exchange  collisions  represented  by  Eqs.  (3.21a,  b,  c) 
need  to  be  retained. 

(  i  ’ 

The  V-V-T  processes  represent  the  most  important  mechanism  for 
determing  the  vibrational  population  distribution  of  the  CO  laser  system. 
Although  excitation  energy  input  comes  from  electron  impact  collisions, 
the  rates  of  V-V-T  resonant  exchange  collisions  are  much  faster,  and 

thus,  the  energy  is  rapidly  redistributed  among  the  vibrational  levels. 

34  35-41 

It  has,  in  fact,  previously  been  shown  by  Treanor,  et  al  and  others 
that  the  kinetics  of  diatomic  molecular  gases  dominated  by  rapid  V-V-T 
exchange  collisions  can  lead  to  highly  excited,  non-Boltzmann  distributions 
of  vibrational  level  populations. 

From  the  principle  of  detailed  balancing,  the  transition  rates  for  each 
of  the  microscopic  reversible  processes  above  will  be  equal.  In  the  case 
of  the  V-T  processes  (3.20)  and  the  V-V-T  processes  (3.21a,  b,  c),  excess 
energy  has  been  transferred  to  (or  from)  the  thermal  reservoir  at  tempera¬ 
ture  T,  and  the  initial  and  final  states  of  the  reservoir  have  been  tacitly 
suppressed  in  writing  these  reactions.  The  heat  bath  can  be  considered 
to  be  a  system  vath  a  large  number  of  degrees  of  freedom,  with  states  |f  > 
that  span  a  continuous  energy  spectrum.  The  assumption  of  thermal 
equilibrium  at  temperature  T  means  that  the  heat  bath  can  be  defined  as 
a  statistical  mixture  with  probabilities  P^  =  exp  (-E^kT)  for  finding  it  in 
any  of  the  states  |  f  >.  (The  assumption  of  infinite  heat  capacity  is  equivalent 
to  saying  that  the  probabilities  P^  do  not  change  if  energy  is  exchanged  with 
the  reservoir;  i.e.,  the  termperature  T  does  not  change.)  If  V-T  processes 
are  written  to  explicitly  include  reference  to  the  initial  and  final  states  of 
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the  reservoir,  then,  for  example 


|  i  >  +  A(r+1)  +  A  2  If  >  +  A(r)  +  A 


(3.22) 


E.  +  E  ,  =  E.  +  E 
1  r+1  f  r 


(3.23) 


If  the  rates  for  the  microscopic  process  (3.  22)  are  denoted  by 

p  .  .  ,  and  p  .  „  .,  the  principle  of  detailed  balancing  requires 

r+1,  i-»r,  f  r,  f-^r+1,  i 

that 


Pr+1,  i-*r,  f  ~  Pr,  f-»r+l,i* 


(3.24) 


However,  since  the  internal  states  of  the  thermal  reservoir  are  not  of 
interest,  the  total  V-T  transition  rates  can  be  obtained  as  follows: 


P 

P 


r+l-«r 

r-*r+l 


=  x5Pi  Pr+1’  Ur‘f 

"  Pf  Pr,  f-*r+l,  i 


(3.25a) 

(3.25b) 


i.  e. ,  the  total  transition  probabilities  are  obtained  by  summing  over  all 
final  lattice  states,  and  over  all  initial  states  weighted  by  the  a  priori  / 
probabilities  p.  for  being  in  states  |i>.  Since  P^  =  P^  exp  £(-E^  +  E^)/kTj  = 
exp  £(Er  -  Er+j,'/kTj,  Eqs,  (3.24)  and  (3.25)  lead  to  the  thermodynamic 
expression  of  the  principle  of  detailed  balancing, 


P 


r+l-»r 


(3.26) 


Similar  equations  hold  for  all  of  the  forward  and  reverse  transition 
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K-*S3W»'5HJ$P^ 


and  fcr  (3.21b)  by 


Qr+1,  R-r,  R+l  =  Qr,  R+l-r+1,  R  *** 


/E  -  E  +  F  -  F 
r  r+l  r  R  R+l 


kT 


)• 


The  rate  of  excitation  of  level  r  due  to  V-T  and  V-V-T  collision 


processes  will  be  given  by 


VT  + 

Rr  (A  ,  A,  =  Pr+Ur 


-  P 


r-»r-l 


VT  + 
Rr 


r-*r-l 


\ 


“r+l  n  *  *XP  \  kT 


/  Er+1  ~  Er\ 

\  kT  /  "r” 

:  /Er-Er-l\ 

nr"  "  CXP  V  ~Tt - /  nr-in 


“r+lN  •  exp  (*  —*f  )nrN 
nrN-cxp(-  '  KTr'1)mr-lN 


(3.27a) 


(3.27b) 


(3.28a) 


(3.28b) 
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RVVT  t  f  .  rp  fn  n  . 

r  r+1,  s-l-»r,  s  r+1  s-l 

L 

t  Er+l  +  Ee-l~ErEs\  1 

exp\“ - kT - /VsJ 


(3.29a) 


Fr,  s-»r-l,  b+1 


T  /  Er  +  E.-Er-l-Es+l)  1 

[“A  -  V - kT - 7  nr-ln8+lJ 

Rm(*,yi.Zo(1,1  B  [■„»,. 

r  R  r+lfR-l-*r,  R  r+1  R-l 

/  Er+1  '  Er  +  FR-1  -  FR  ^  „  1 

“p  V  - kT - KNR  J 


-  7.  x 

Pk  r,  R-*r-l,  R+1 
[nrNR  '  exp  ( 


(3.29b) 


E  -  E  +  F  -  F 
r  r-l  R  R+1 


)n  I 

r-l  R+1  J 


Before  discussing  the  rates  in  detail,  it  is  instructive  at  this  point  to 

illustrate  the  physical  significance  of  fast  V-V-T  cross -relaxation  with 

42 

a  simple  example.  Javan  has  shown  how  the  interplay  between  two 
nearly  resonant  2-level  systems  in  contact  with  thermal  reservoirs  at 
different'  temperatures  (e.  g. ,  the  electron  and  translational  baths)  can 
lead  to  interesting  equilibrium  population  distributions  under  different 
conditions,  depending  upon  the  relative  magnitudes  of  the  rates  of  pumping 
and  relaxation  mechanisms.  Consider  the  two  systems  shown  in  Fig.  3.  3; 
system  s  has  populations  n^  and  n^  in  levels  a  and  b,  which  are  separated 
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System  s  System  S 


Figure  3.3a.  Two  simple  2-level  systems  interacting  with 
resonant  exchange  collisions. 


by  an  energy  e,  while  system  S  has  populations  and  in  levels  1  and 
2,  which  are  separated  by  an  energy  E.  Assume  that  level  b  is  pumped 
at  a  rate  r,  and  level  2  with  a  rate  R,  and  that  near  resonant  exchange 
collisions  between  the  two  systems  are  possible.  If  the  translational 
mode  is  in  thermodynamic  equilibrium  at  temperature  T,  the  collisional 
exchange  rates  will  be  related  by 


P,,2-.b.l  *  <4E/kT)  Pb.l-».2  >  <3’30> 

where  AE  =(E  -  e),  and  the  rate  equations  for  the  populations  n^  and 
will  be  given  by 


dn^/dt  =  r  +  C 
dN2/dt  =  R  -  C, 


(3.31) 


where  the  collision  term  C  is 


C 


n  N 

a 


2P», 


2-b,l 


(3.32) 


If  it  is  assumed  that  the  coll;  ion  term  C  is  the  dominant  rate  for  system 
s,  then  r  can  be  neglected,  *nd  in  the  steady  state. 


<VN1>  <Pa,2.b,l/Pb,l.a, 


2)  =  (N2/N1)exp  (AE/kT). 


(3.33) 


On  the  other  hand,  suppose  that  the  pumping  term  R  is  the  dominant  rate 
for  system  S,  with  collisional  loss  terms  negligible,  and  assume  that  R 
brings  the  populations  N  to  a  temperature  ®: 


N2/Ni  =  exp  (  -  E/k©). 


(3.34) 
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For  example,  system  S  could  be  thought  of  as  being  pumped  by  electrons 

to  a  high  temperature,  ©  =  T  .  Eq.  (3.  33)  for  the  populations  of  system 

e 

s  becomes 


n^/na  =  exp  (AE/kT  -  E/kT^) 

Eq.  (3.35)  has  a  very  interesting  structure,  for  notice  that,  if 


(3.35) 


(E  -  e)/T  >.E/T  , 
e 


(3.36) 


then  n^/n^>  1,  corresponding  to  a  situation  with  total  inversion.  That 
is^  because  of  the  energy  difference  AE  =  (E-e),  system  S  can  pump 
system  s  to  a  totally  inverted  condition  in  the  situation  where  vibra¬ 
tional  cross -relaxation  is  fast  compared  with  the  other  pumping  and 
relaxation  rates  of  system  s.  If  a  temperature  $  is  defined  for  system 
s,  it  follows  that 


e  E  E-e 

B  ©  =  T 


(3.37) 


A  plot  of  this  equation  is  shown  in  Fig.  3.  3,  which  shows  that  for  values 
of  ©which  are  sufficiently  large,  the  corresponding  temperature  B 
becomes  negative.  Clearly,  the  effect  becomes  more  pronounced  as 
the  kinetic  temperature  T  is  decreased,  since  the  Boltzmann  factor 
exp(4E/kT)  that  expresses  the  ratio  between  the  forward  and  reverse 
processes  increases.  That  is,  for  lower  translational  temperatures, 
system  S  can  pump  system  s  more  efficiently.  Schematically,  this  can 
be  seen  in  Fig.  3.3,  in  which  the  vertical  asymptote  line  would  shift 
toward  the  left  as  T  is  decreased. 
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Based  upon  this  simple  example,  it  is  apparent  how  rapid  V-V-T  cross¬ 
relaxation  can  have  such  profound  significance  for  the  kinetics  of  the 
vibrational  populations  in  the  CO  laser  system.  Since  has  larger 
energy  spacings  than  CO,  can  pump  the  CO  levels;  furthermore, 
since  the  upper  vibrational  level  spacings  of  CO  decrease  due  to 
anharmonic  effects,  molecules  in  the  lower  levels  of  CO  can  pump  mole¬ 
cules  which  are  in  the  upper  levels  to  higher  vibrational  temperatures. 
This  is  the  physical  origin  of  the  so-called  "anharmonic  pumping" 
mechanism. 

43-66 

There  are  several  theories  of  V-T  and  V-V-T  collision  rates,  and 

thus,  the  dependence  of  these  transition  probabilities  on  vibrational  level 

and  translational  temperature  is  one  of  the  greatest  uncertainties  for 

constructing  a  suitable  molecular  kinetic  model.  The  form  of  the  rates 

43-45  ... 

to  be  used  he  re  is  taken  from  a  modified  form  of  SSH  theory  with  tne 
V-V-T  transition  rates  given  by 

2* 2 1  T  €  r 

Pr,s-l-r-l,s  =  0 J  h°  eXp(T) 

exp  (  0rs/2T)  F(gj 

where  x^  was  the  anharmonicity  factor  defined  by 
kinetic  collision  frequency,  6  =  (E  -  E  ,  -  E 

a  parameter  originating  from  the  use  of  a  Lenard 
potential.  The  quantities  b  and  g^  are  defined  b} 

b  =  lb*  4 Mi  2  k/h2 


(nfc-) 
) 


(3.38) 


rs 


Eq.  (2.22),  f  is  the 
c 

+  E  ,)/k,  and  €  is 
i  s-l 

-Jones  interaction 


r 


(3.39) 
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and 

gX*b0  2/8T  (3.40) 

where  Mis  the  reduced  mass,  and  /  is  a  measure  of  the  effective  inter¬ 
action  length,  and  is  typically  ~.  2  X,  which  we  shall  assume  here.  If  fl 

o  2 

is  expressed  in  molecular  weight,  and  t  is  given  in  A,  then  b  =  .  8153  Ml  . 

The  function  F(g  )  is  an  expression  of  the  "adiabaticity"  of  the  collision, 
r  8 

00 

F(g)  =  g 2  J  d*e"*csch2  (g/y/T ),  (3.41) 

and  which  becomes  approximately 

F(g)  *  8(>r/3)1/2  g?/3  exp  (  -  3g2/3)  (3.42) 

in  the  limit  that  g  is  large  (i.e. ,  off  resonant).  A  modification  of  the 

usual  SSH  theory  (which  uses  Eq.  (3.42)  for  all  cases'  was  suggested 

59  58 

by  Keck  and  Carrier  and  Bray  who  use  an  expression 

F(g)  =  -j  exp  (  -  2g/3)  (3  -  exp  (-  2g/3)]  ,  (3.43) 

which  is  approximately  valid  (to  within  about  20%)  for  all  values  of  g<  20. 

For  the  calculations  to  be  performed  here,  it  is  the  structure  of  Eq.  (3.38) 
which  is  important,  since  it  is  necessary  to  have  some  assumptions  for 
calculating  the  V-V-T  rates  for  high  vibrational  levels  in  terms  of  the 
rate  which  shall  be  assumed  to  be  known.  Thus,  Eq.  (3.38) 
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can  be  written  as 


P  =  P  rs 

r,  s-l-»r-l,  s  10-01 


exp  (  Q^JZT)  -  |3  -  exp  (  -  2grs/3)]  exP  (  -  2grB/3)  (3-  44) 

Likewise,  the  rates  for  V-T  processes  can  be  taken  from  SSH  theory, 
and  are  given  by 

pr-r.i =  p-o r  (nr-r)  «p  <  ,'2T>  *•«*,»/«*,)  <3- 45 » 


where  0  = 
r 


=  (E  -  E  .  )/k,  $  =  $  -6  ,  and  g  =  b$  1 

r  r-1  r,l  r  1’  r  r 


/8T. 


Experimental  rates  for  V-T  processes  have  been  measured  by  several 
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workers,  and  results  have  been  summarized  by  Millikan  ;  some  of 

his  results  are  shown  in  Figure  3.4  and  3.5  .  V-V-T  rates,  on  the 

other  hand,  are  typically  somewhat  scarce,  and  those  that  are  reported 

are  usually  measured  at  high  temperatures  (Z  1000°K),  characteristic 

64 

of  shock  tube  experiments. 
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Hancock  and  Smith  and  Yardley  have  made  some  experimental  meas¬ 
urements  of  CO-CO  vibrational  exchange  rates,  and  we  are  presently 

making  preparations  to  infer  some  of  these  rates  from  double-probe  and 

21 

double -Q  switch  experiments,  to  be  described  in  a  companion  report 

SSH  theory  employs  the  assumption  that  Bhort-range  repulsive  interactions 

(e.g. ,  the  Lenard-JoneB  potential)  are  the  dominant  mechanism  for  V-V-T 

46 

energy  exchange.  However,  as  Mahan  has  suggested,  molecules  with 
large  transition  dipole  matrix  elements  may  be  able  to  transfer  vibrational 
quanta  more  efficiently  by  means  of  resonant  processes  which  result  from 
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49-53 

long  range  dipole-dipole  interactions.  Sharma  has  made  several 

calculations  of  the  energy  transition  rates  for  dipole  and  multiple  coupling, 

and  has  shown  that  the  temperature  dependence  of  the  transitions  may 

differ  markedly  from  the  results  that  SSH  theory  predicts.  According  to 

66 

Hancock  and  Smith  ,  calculations  for  CO-CO  based  on  a  model  similar 

to  Sharma *s  have  been  performed  by  Morley  and  Smith,  and  will  be  published 

67 

later.  Jeffers  and  Kelley  have  also  recently  submitted  results  of  calculations 
based  on  both  long  range  and  short  range  interactions.  As  more  reliable 
data  becomes  available,  it  will,  of  course,  be  possible  to  make  improvements 
in  the  kinetic  model  for  CO.  At  the  present  time,  the  SSH  theory  is  employed, 
with  tfce  calculations  to  be  presented  later  to  be  based  on  Eq.  (3.44)  and  (3.45) 
for  reasonable  values  of  etc* 

3.  5  General  Discussion.  The  vibrational  levels  of  A  and  B  are  strongly 
coupled  among  themselves,  and  to  each  other,  through  the  V-V-T  exchange 
collisions.  The  situation  that  will  result  for  the  population  distribution 
in  the  steady  state  will  depend,  as  the  master  equation  (3.16)  indicates, 
upon  the  relative  importance  of  many  different  pumping,  relaxation,  and 
cross -relaxation  rates.  In  this  section,  a  somewhat  qualitative  discussion 
of  the  molecular  kinetics  will  be  given,  with  the  results  of  computer 
calculations  to  be  presented  subsequently.  In  spirit,  the  discussion  here 
will  follow  the  same  approach  adopted  above  for  the  simple  model  of  the 
two  interacting  2-level  systems:  viz, ,  if  any  term  R^  in  Eq.  (3.1b)  for 
two  coupled  diatomic  species  A  and  B  is  dominant  (with  all  remaining  terms 

X  X 

negligible  relative  to  Rf),  then  =  0  will  be  the  most  important  condition 
for  determining  the  form  of  the  equilibrium  population  distribution. 

a)  V-T  Collisions.  If  V-T- collisions  are  dominant,  the  largest  term  in 

VT  ♦  VT  f 

Eq.  (3.1b)  will  be  R  (A  ,  A)  or  R  (A  ,B)  or  both.  If  all  other  rates  are 

r  r 

neglected  in  comparison,  then  the  equilibrium  vibrational  distribution  will 
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be  defined  by  the  set  of  equations 


/friLlfr.) 

'  kT  ' 

(nP  A1  +  NQ  )  (n  -  e  A  n  ) 

r+l-»r  r+l-*r  r+1  r 


=  (nP  .+  NQ  .)  (n  -  e 
r-*r-l  r-*r-l  r 


,E  -  E  ,  v 

£ - Li) 

'  kT  / 


(3.46) 


n 

r-1 


For  r  =  0,  only  one  term  enters,  viz., 


kT 


<nPi~o ♦  NCW  <"1  • e  V  ■  °- 


(3.  47) 


Thus,  the  vibrational  distribution  becomes,  under  this  limiting  condition, 


n 

r 


=  n  e 
o 


-Er/kT 


(3.48) 


1.  e. ,  the  vibrational  levels  achieve  thermal  equilibrium  at  the  wall 
temperature  T.  Since  V-T  rates  increase  with  quantum  level  (both 
because  of  the  larger  transition  matrix  element,  and  the  smaller  amount 
of  energy  4E  that  is  exchanged  with  translation),  these  processes  have  a 
cooling  effect  upon  the  higher  levels. 


b)  Electron  Collisions.  If  the  electron  pumping  rate  Re,  given  by 
Eq.  (3.18),  is  the  dominant  term  in  Eq.  (3.1b),  then  the  population 
distribution  will  be  determined  mainly  by  setting  =  0,  which  leads  to 

n  =  n  exp  (-E  /kT  ) . 
r  o  r  r  e 

That  is,  if  electron  collision  rates  are  such  that  other  relaxation  and 
cross -relaxation  mechanisms  are  negligible  in  comparison,  then  the 
vibrational  levels  will  be  raised  to  a  temperature 
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c)  V-V-T  Processes.  In  the  case  of  V-V-T  processes,  each  set  of 

levels  (of  both  systems  A,  B)  is  in  a  situation  where  it  attempts  to 

equilibrate  with  itself,  and  each  set  is  in  communication  with  the  other. 

t  t«  t  t 

If  V-V-T  processes  of  the  type  A  +  A*-  A  +  A  dominate  the  steady  state 
for  species  A,  and  all  other  mechanisms  (including  A*  +  A*+  B*  )  are 

negligible,  then  the  equilibrium  state  for  the  vibrational  levels  of  A  is 
defined  by  IT  (AT,  AT  )  =  0: 


^  Pr+1.  s-l-r,  s, 


(E  —  E  +  E  ,  —  E  \ 
r+1  r  s-l  si 

■ - — - Jn  n 

kT  /  r  s 


n  n 
r  s 


“  ^  Pr*  8 

a» 

/  Er-Er-l  +  E.“E.+l\ 

CXP  [- - kT - /  nr-l“ 


(3.49) 


It  can  be  verified  by  direct  substitution  that  a  solution  to  this  equation 
is  given  by 


n  =  n  e“r y  exp  (-E  /kT), 
r  o  r 


(3.  50) 


where  y  is  a  single  parameter  for  all  of  the  energy  levels  r(r  =  0,  1,  2. . . ), 
Likewise,  if  V-V-T  processes  for  species  B  dominate  all  other  (including 
even  V-V-T  collisions  of  the  type  B^  +  /$?Bt+  it ),  then  the  distribution  of 


populations  will  be 


\ 


NR=Noe'Rrexp(-FR/kT). 


(3.  51) 


In  general,  y  need  not  equal  /7  since  the  two  systems  equilibrate 
independently.  If  vibrational  temperatures  9^  for  each  level  r  of 
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species  A  are  defined  by 


n  =  n  exp  (-  E  /k0  ), 
r  o  r  r  * 


r  o 

r  r 

then 

E 

E 

«• 

A 

ii 

kT  +  ry- 

For  r  =  1, 

we  get 

r_5_ 

E1 

“I 

kT  * 

so  that 

9t  _ 

frY  ( 

'I!l  .A l-1. 

T 

k«.  ' 

^Er  /J 

The  distribution  function  (3.  50)  becomes 

/  rEl\  /rEl-Er\ 
nr  =  noCXp  k^j  CXP  V""kT“7* 


(3.52) 


(3.53) 


(3.54) 


(3.55) 


(3.  56) 


Clearly,  if  the  levels  r  describe  harmonic  oscillators,  then  Ef  =  rE^, 
and  the  distribution  of  vibrational  levels  will  be  Boltzmann,  with  some 
temperature  0,  that  characterizes  the  entire  set  of  levels.  Likewise  for 
species  B.  We  shall  assume  for  the  moment  that  each  species  A,  B  is 
described  by  a  harmonic  oscillator,  since  it  will  then  be  possible  to 
discuss  the  effects  of  V-V-T  dominance  in  terms  of  the  convenient 
qualitative  concept  of  a  single  vibrational  temperature  for  each  species. 
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Then 


and 


E1 

y  -  — - 

7 

Ei , 
kT  ’ 

(3.57a) 

„  F1 

Fi 

r  =  _L 

k 

1 

kT  ’ 

(3.57b) 

n  =  n  exp  (-rE,  /k0), 
r  o  1 

(3.  58a) 

N  =  N 

R  o 

exp  (-RFj/k©  ) . 

(3.58b) 

The  result  (3.  57)  was  obtained  by  assuming  that  V-V-T  collisions  between 
A  and  B  of  the  type  (3.  21b)  are  negligible  compared  to  V-V-T  collisions  of 
the  type  (3.21a,  c).  This  would  typically  be  the  case  if  the  levels  of  species 
B  do  not  have  a  good  resonance  match  with  those  of  species  A,  since  V-V-T 
collision  rates  drop  rapidly  with  an  increase  in  the  energy  defect  that  must 
be  exchanged  with  the  translational  mode.  In  general,  therefore,  the 
vibrational  temperature  6  and  ©attained  by  the  levels  of  A  and  B  under  a 
dominance  of  V-V-T  collisions  of  types  (3.  21a,  c)  would  not  be  equal. 

Alternatively,  if  we  assume  that  all  V-V-T  processes  have  rates  of  the 
same  order  of  magnitude,  then  the  equilibrium  distribution  attained  by 
the  set  of  levels  of  A  in  a  steady  state  dominated  by  exchange  collisions 
will  be  determined  instead  from 

RyVV,  Af  )  +  RyVT(A* ,  Bf  )  =  0  (3.  59) 
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..-spajf, 


It  can  be  shown  that,  under  this  circumstance,  the  solution  for  the  level 

populations  n  and  N_  of  the  two  species  is  given  by 
r  R 


n 

r 


(3.60a) 

(3.60b) 


where  y  is  the  same  for  each  set  of  levels.  Thus  for  harmonic  oscillators, 
if  all  the  V-V-T  rates  are  comparable,  the  solutions  (3.58)  give  the  popu¬ 
lation  distributions  subject  to  the  constraint 


E 


1 


e 


(3.61) 


This  equation  is  obtained  from  (3.  57)  by  setting  y  =  /"*,  and  is  identical 
to  Eq.  (3.37),  which  was  derived  earlier  for  the  simple  two-level  model. 


In  paragraphs  a),  b),  and  c)  above,  it  has  been  shown  how  the  vibrational 
levels  can  achieve  quire  different  distributions,  depending  upon  which  of 
the  relative  rates  for  several  competing  mechanisms  is  dominant.  In 
actual  practice,  the  equilibrium  distribution  is  determined  by  an  interplay 
of  these  processes,  and  typically,  vibrational  temperatures  of  5.10,  000°K 
are  attained.  For  the  higher  quantum  levels,  say  r  .>25,  V-T  processes 
dominate  V-V-T  processes,  and  the  upper  levels  can  be  expected  to  reach 
lower  local  temperatures  by  increased  contact  with  the  wall  reservoir. 


It  was  pointed  out  previously  that,  due  to  the  smaller  energy  spacing  between 
the  upper  vibrational  levels  of  CO,  the  fast  V-V-T  cross  relaxation  processes 
can  excite  the  upper  levels  to  even  higher  vibrational  temperatures  by  the 
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process  of  "anharmonic  pumping.  "  N2  has  levels  which  are  slightly  greater 

than  those  for  CO,  and  can,  therefore,  pump  the  CO  levels.  Since  N_  does 

£ 

not  have  a  very  fast  V-T  relaxation  rate,  it  can  serve  as  an  efficient 
reservoir  for  storing  vibrational  energy,  which  is  subsequently  transferred 
to  CO.  (Addition  of  N2  to  an  electric  discharge  system  also  has  important 
consequences  for  the  plasma  characteristics,  which  will  be  discussed  later.) 


Further  interesting  conclusions  can  also  be  obtained  by  considering  rates 
other  than  the  dominant  ones  for  the  system.  For  example,  it  is  possible 
to  illustrate  the  effects  of  N.,  and  CO  V- V  pumping  by  showing  how  a 
"vibrationally  cold"  system  can  transfer  energy  to  a  "vibrationally  hot" 
one,  which  may  seem  paradoxical  thermodynamically.  Consider  the  situation 
where  the  two  species  A,  B  attain  quasi-Boltzmann  equilibrium  distributions 
separately,  with  the  V-V-T  exchange  collisions  between  A  and  B  negligible. 

(I.  e. ,  suppose  that,  perhaps  because  of  severe  energy  defect,  the  off- 
resonance  rates  for  cross -relaxation  between  systems  A  and  B  are  much 
slower. )  Then  on  a  short  time  scale,  the  vibrational  populations  attained  by 
the  two  species  will  be  expressed  by  (3.  50)  and  (3.  51),  where  y  and  Tare  not, 
in  general,  equal.  On  a  longer  time  scale,  the  two  systems  will  interact, 
with  the  rate  of  energy  change  in  species  A  due  to  V-V-T  exchange  with  B 
given  by 


Z  ERVVV.B»>  * 


r  r 


R,r 


r+1,  R-l-*r ,  R 


n  ,,N_ 
r+1  R 


’  ]£rErQr,R-r-l,R+l 


_j  -  exp 
nrNR  '  exP  ( 


(3.62) 


E  -  E  +F 
r+1  r 


kT 


E  -  E  +  F  -  F 
r  r-1  R  R+1 
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-)n  ,N_ 
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Assume,  for  simplicity,  that  the  rates  are  determined  by  a  dipole - 

dipole  interaction  between  harmonic  oscillators,  and  are  given  by 

Q  «  .  — r(R+l)Q,  _  . 

r,  R-»r-l,  R+l  1,  0-*0, 1 


Eq.  (3.62)  becomes,  upon  substituting  the  quasi>steady  distributions 
(3.50),  (3.51) 


VVT  t  t  Fj/kT  V"* 

rRr  (*.B,)-E1Q1^(W,1.  er,  ^,RNrL 


rn 


(3.63) 


If  this  expression  is  negatr  e,  then  energy  is  being  transferred  from 
A  to  B  through  the  V-V-T  exchange  mechanism  of  Eq.  (3.21b).  This 
condition  will,  therefore,  occur  when /*>y,  i.  e. ,  for 


T 


(energy  flow:  A 


B) 


(3.64) 


Alternatively,  if  in<y,  the  expression  (3.63)  is  positive,  indicating  that 
energy  is  being  supplied  to  system  A  by  system  B: 


(energy  flow:  B 


A) 


(3.65) 


Figure  3.  6  shows  a  schematic  plot  of  the  hyperbola  defined  by  (Ej/0)  - 
(Fj/©)  *  (Ej  -  Fj)/T  for  the  case  E^<  F^.  In  the  region  above  this  curve, 
energy  exchange  is  from  A  to  B,  while  in  the  region  below,  it  is  from 
B  to  A.  The  straight  line  corresponds  to  equal  vibrational  temperatures, 
0  =©.  Note  that,  in  the  shaded  region  of  the  curve  0>®but  the  energy 
flow  is  from  B  -*  A.  I.  e. ,  under  this  unique  situation,  energy  can  flow 
from  the  "colder"  system  to  the  "hotter"  system. 
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0  (Vibrational  Temp,  of  A) 


0  (Vibrational  Temp,  of  B) 


Figure  3.  6.  V-V-T  cross  relaxation  for  two  harmonic 

diatomic  species  A  and  B,  with  energy  spacings 
Ej  and  Fj(  and  vibrational  temperatures  0  and®, 
respectively.  The  shaded  region  corresponds 
to  the  case  0  >  ®,  and  thus  represents  a  situa¬ 
tion  where  the  flow  of  energy  (B*A)  is  from  the 
"colder"  to  the  "hotter"  system.  (Adapted  from 
Reference  35). 
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3.6  Computer  Calculations  and  Numerical  Results.  The  computer  program 
that  has  been  developed  to  analyze  the  steady  state  is  capable  of  accepting 
all  of  the  above  inputs,  for  arbitrary  pressures,  gas  mixtures,  temperatures, 
etc.  At  the  present  time,  the  populations  of  a  second  diatomic  species  such 
as  N2  (if  present)  are  treated  parametrically  as  specified  input,  rather  than  as 
independent  variables  to  be  included  in  a  second  coupled  set  of  master  equa¬ 
tions.  (However,  if  and  when  it  becomes  necessary,  the  program  can  be 
easily  extended  to  treat  the  populations  of  both  of  two  interacting  diatomic 
species  on  an  equal  basis. )  The  calculation  can  be  carried  out  for  an 
arbitrary  number  of  CO  levels  (in  practice,  usually  £40  is  sufficient),  and 
will  produce  all  of  the  relevant  information  needed  for  predicting  the 
characteristics  of  the  laser  amplifier  medium.  The  output  produced  consists 
of  the  vibrational  population  and  temperature  distribution,  the  rates  of  energy 
transfer  through  the  competing  mechanisms  present,  laser  gains,  saturation 
intensities,  stimulated  energy  extraction,  and  quantum  efficiency.  Thus,  it 
is  possible  to  obtain  small-signal  or  saturated  gains  by  appropriate  input 
intensities  for  either  single  or  multiline  operation. 

3.6. 1  Parametric  Calculations.  In  order  to  perform  a  preliminary  check  on 
the  computer  program,  some  initial  special  cases  for  which  the  solution  is 
known  in  advance  were  run  as  part  of  the  debugging  procedure.  As  was  pointed 
out  in  Section  3.5,  turning  off  all  processes  except  for  the  electron  pumping 
should  result  in  a  distribution  for  which  all  of  the  vibrational  temperatures 
are  equal  to  Te,  and  this  result  was  correctly  verified.  Alternatively,  by 
turning  off  everything  but  V-T  and  V-V-T  relaxation,  the  system  should 
approach  a  steady  state  for  which  all  vibrational  temperatures  equilibrate 
to  the  wall  temperature  T,  and  this  was  also  confirmed.  By  turning  off  the 
anharmonicity  of  the  diatomic  species,  an  equilibrium  distribution  whose 
vibrational  temperatures  satisfy  relation  (3.  61)  should  result,  providing 
other  processes  are  negligible,  and  this  also  was  correctly  obtained. 
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To  assess  the  effects  of  a  variation  of  parameters  on  the  resulting  popula¬ 
tion  distribution  in  the  CO  laser,  a  variety  of  comparative  calculations  were 
performed  for  a  low  pressure  system,  using  reasonable  assumptions  for  the 
relaxation  rates  and  values  for  electron  pumping  rates  that  would  be  typical 
of  an  electrically  excited  longitudinal  discharge  configuration.  In  order  to 
construct  a  reliable  model,  it  will  be  necessary  to  perform  measurements 
of  the  small  signal  gain  coefficient  for  several  controlled  sets  of  operating 
parameters,  and  experimental  work  toward  this  objective  is  now  in  progress. 

In  Figure  3.7  ,  the  effects  of  increasing  partial  pressure  of  CO  illustrate 

how  anharmonic  pumping  can  cause  the  excitation  to  higher  v  levels.  From 
a  different  point  of  view,  Figure  3.8  shows  how  a  variation  in  the  rate 
of  V-V-T  cross -relaxation  can  have  a  profound  effect  on  the  resulting 
vibrational  distribution.  In  Figure  3.9  the  interaction  length  t ,  which 
is  critical  parameter  which  measures  how  fast  the  off -resonance  V-V-T 
rates  decrease,  has  been  varied.  Figure  3.10  indicates  the  effect  of 
increasing  the  electron  number  density  to  increase  the  pumping  rate. 

Figure  3. 11  corresponds  to  addition  of  N£»  whose  vibrational  temperature 
is  parametrically  fixed  at  two  different  constant  temperatures. 


3.6.  *  Energy  Transfer  and  Quantum  Efficiency.  Calculation  of  the  rates 
of  energy  transfer  for  all  of  the  kinetic  and  radiative  mechanisms  pre¬ 
sent  in  the  CO  laser  amplifier  medium  are  important  for  predicting  detailed 
rates  of  heating,  electrical  energy  input,  single  and  multiline  extraction 
capacity,  and  quantum  efficiency.  The  rates  of  V-T  and  V-V-T  heating  are 
important  parameters  for  design  considerations  of  high  power  CO  laser 
devices.  The  quantum  efficiency,  which  can  be  defined  as  the  percentage 
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Figure  3.  9.  Effect  of  Variation  of  Interaction  Length  l  of  SSH  theory 
on  CO  Vibrational  Population  Distribution. 
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Figure  3. 11.  Effect  of  N2  (included  parametrically  at  fixed  vibrational 
temperature)  on  CO  Population  Distribution. 


81 


of  coherent  radiation  extracted  for  a  given  electrical  input,  is  a  function 
of  the  input  line  intensities  and  can  be  enhanced  (approaching  the  theoretical 
limit  of  100%)  if  sufficient  numbers  of  lines  with  intensities  near  and  above 
saturation  are  included.  Saturation  intensities  for  single  or  multiline 
operation  can  be  calculated,  as  well  as  volumetric  extraction  energies, 
for  each  line  present.  The  steady-state  rates  of  energy  transfer  for  each 
of  the  process  considered  in  Eq.  (2. 16)  are  summarized  below: 


a.  Electrical  Energy  Input.  The  rate  of  electrical  excitation  is  given  by 

U*=  £  E  [n  R*  -nR*  ]  (3.66) 

r  r  r-*s  s  8«#r 

r,  s 

b.  V-V-T  Heating.  For  the  process  CO(r)  4  CO(s-l)  -»CO(r-l)  4-  CO(s)  4 

4e  ,  where  an  amount  of  energy  AE  =  (E  -E  J-(E  -E  )  has  been 
r ,  s  r , s  rr— lss— l 

transferred  to  kinetic  energy,  the  rate  of  V-V-T  heating  due  to  the  anharmonic 
defect  is  given  by 


U 


VV 


S(E  -E  ,)[n  n  ,P  .  .  -  n  ,n  P  .  .]  (3.67) 

r  r-1  r  s-1  r,s-l-»r-l,  s  r-1  s  r-1,  s-»r,s-l 


If  there  is  another  diatomic  species  B  in  the  gas  mixture,  then  the  rate  of  heat 
generation  for  the  process  CO(r)  4-  B(R-l)-*  CO(r-l)  4-  B(R)  4-  ^E^  ^  where 

*Er.R  =  <Er-Er.l>  '  <FR-FR.1>-  “  given  b'r 


UV  ’  D  <Er-Er-l-FR  +  FR-l)  ["rNR-1  Qr  R., 
r,R 


"nr -1  NR  Qr-l,R-*r,  R-l] 


(3.68) 
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c. 


V-T  Heating.  For  V -T  collisional  relaxation  processes  such  as 

CO(r)  +  X  -►CO(r-l)  +  X  +  AE  ,  where  A E  =  E  -E  ,  the  rate  of 

r  r  r  r-1 

V-T  heating  is  given  by 


U 


VT 


.  £ 


(E  -E  .)  [n  P  .  -  n  .  P  .  ]  N 

'  r  r-1  r  r*^  r-1  r-1  r-1-*  r  3 


(3.69) 


d.  Vibrational  Energy  Transfer.  For  a  situation  where  there  is  a  second 
diatomic  gas  B  included  in  the  mixture,  there  will  be  a  transfer  of  enorgy 
from  one  vibrational  species  to  the  other.  The  rate  of  energy  flow  from  B 
to  CO  by  the  process  CO(r)  +  B(R-l)  -*  CO(r-l)  +  B(R)  +  4  E  is  given  by 


uB"CO*-rVFR-F 


_  .)  [n  N_  .Q  _  .  ,  _ 

R-1  r  R-1  r, R-1 r-1,  R 


*nr-lNRQr-l,  R-*r,  R-l] 


(3.70) 


or  equivalently,  the  transfer  of  vibrational  energy  from  CO  to  B  is 

*CO-*B  *B-CO  ;TCO-*BiT'TVV 

U  =  -  U  .  It  is  easy  to  verify  that  U  rePre8en;s 

the  total  loss  of  vibrational  energy  from  CO  in  a  V  -V-T  exchange  collision 

with  B. 


e.  Spontaneous  Radiation.  Spontaneous  radiation  processes  represent 
a  pure  loss  mechanism,  and  thus,  their  presence  limits  the  total  attainable 
quantum  efficiency  to  a  value  less  than  100%.  For  the  Ay  =  1  and  2  decay 
processes  considered  in  the  present  model,  the  rate  of  radiative  energy  loss 

is 


U 


■P 


=  ^  n  (E  - 


E  .)  A  ,+ 
r-1  r-+r-l 


Sn  (E  - 


Er-2)  Ar-*r-Z 


(3.71) 


f.  Stimulated  Emission.  If  stimulated  emission  processes  are  important, 
i.  e. ,  if  the  intensities  of  the  input  lines  are  comparable  to  or  greater  than 
the  saturation  parameters,  the  gains  for  these  transitions  will  saturate. 
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The  extraction  power  per  unit  volume  will  be  given  by 


dl/dx  =  a(I)I 


(3.72) 


where  I  is  the  input  intensity,  and  a(l)  is  the  saturated  gam.  The  total 
stimulated  extraction  from  the  amplifier  medium  will  be  a  sum  of  terms 
(3.  72)  for  all  of  the  input  laser  lines  present.  As  an  alternative  calculation, 
the  total  rate  of  energy  transferred  to  coherent  amplified  radiation  would  be 
given  by 


®  ,tto>  ■  ?  <Er  '  Er-l>[\Sr  ^  r-1  ‘  Vl  Sr-WrJ  <3'  73> 


From  an  initial  calculation  of  the  unsaturated  amplifier  characteristics, 
for  which  no  stimulated  emission  processes  are  included,  the  small-signal 
gains  can  be  obtained.  Saturation  intensities  can  then  be  calculated  for 
single  or  multiline  operation  from  an  asymptotic. limit, 


I  =  lim  a(  I  )  I /a  (3.74) 

s  _  o 

I«*oo 

where  a  (I)  and  a  are  saturated  and  small -signal  gains  for  any  desired 
o 

transition.  In  practice,  this  involves  merely  that  an  input  value  be  used  for 

8  2 

I  which  is  very  large  (e.  g. ,  10  watt/cm  ).  The  complete  generality  available 
for  the  computer  code  input  description  makes  it  possible  to  calculate  many 
interesting  performance  characteristics  of  the  amplifier  medium  --  e.  g. , 
the  cross -saturation  effects  that  determine  the  energy  available  for 

3 

extraction /cm  on  one  transition  in  a  situation  where  other  lines  are 
simultaneously  being  driven  to  saturation.  Knowledge  of  this  sort  is 
extremely  important  for  extrapolating  design  considerations  to  a  high 
energy  CO  laser  system,  where  it  is  imperative  to  extract  as  much  useful 
energy  as  possible  from  several  simultaneous  transitions  with  desirable 

i: 

characteristics. 
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In  the  steady  state,  the  rate  of  electrical  energy  input  into  the  vibrational 
levels  must  be  balanced  by  all  of  the  energy  loss  mechanisms  which  are 
present.  For  the  case  of  an  unsaturated  amplifier,  for  which  all  stimulated 
emission  processes  are  neglected,  the  dominant  mechanism  for  dissipating 
the  input  power  is  by  means  of  V-V-T  heating  which  results  from  the 
enharmonic  defect.  Since  the  V-V-T  rates  as  well  as  the  anharmonic 
defect  increase  for  larger  v  levels,  the  population  distribution  will  be 
flumped  up  far  enough  that  heating  rates  can  balance  electrical  excitation 
rates.  Thus,  if  the  electrical  pumping  is  increased,  the  population 
distribution  will  be  extended  further  up  the  vibrational  ladder.  When 
external  lines  are  introduced  into  the  amplifier,  with  intensities  comparable 
to  or  greater  than  saturation,  stimulated  emission  processes  enter  into 
the  steady  state  energy  balance  condition,  so  that  excitation  energy  is 
dissipated  both  by  radiative  extraction  and  V-V-T  heating  processes.  In 
order  to  maintain  the  population  distribution,  of  course,  V-V-T  heating  is 
always  present,  since  the  vibrational  exchange  collisions  are  the  dominant 
process  for  the  molecular  kinetics.  Thus,  V-V-T  heating  effects  limit  the 
attainable  quantum  efficiency  of  conversion  of  electrical  energy  to  coherent 
radiation,  and  can  also  present  difficult  heating  problems  for  high  power 
devices.  By  operating  the  amplifier  for  low  v-bands,  the  quantum  efficiency 
is  best,  and  it  can  be  improved  by  saturated  extraction  on  several  v-bands, 
Preliminary  calculations  indicate  the  possibility  of  achieving  quantum 
efficiencies  of  70  -  90%  for  high  pressure  CO  laser  devices. 

It  has  been  observed  theoretically  that  the  introduction  of  several  P(  J) 
lines  corresponding  to  a  single  v-band  into  a  CO  laser  amplifier  (assuming 
that  the  assumption  of  rapid  rotational  cross  relaxation  is  valid)  can  result 
in  interesting  energy  transfer  possibilities  when  the  incident  line  intensities 
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are  comparable  to  or  above  the  saturation  intensity.  Exploitation  of  this 
fact  may  provide  the  possibility  for  rotational  line  selection,  since  it  is 
possible  to  adjust  the  characteristics  of  the  amplifier  medium  in  such  a  way 
that  energy  from  one  rotational  line  is  completely  transferred  to  another. 

The  line  selection  device  could  be  operated  at  low  enough  pressures  to  insure 
that  the  lines  provided  by  an  oscillator  enter  with  intensities  above  saturation. 
Indications  of  the  calculations  are  that  the  transfer  can  be  relatively  efficient, 
and  that  heating  difficulties  are  relatively  modest,  which  would  appear  to  be 
advantageous  for  practical  control  of  such  parameters  as  temperature. 
Basically,  this  is  analogous  to  a  parametric  effect,  in  which  the  CO  medium 
is  pumped  by  one  resonant  transition  with  resulting  enhanced  characteristics 
for  another. 


3. 6„  3  Numerical  Technique.  The  system  of  molecular  kinetic  equations 
for  the  steady  state  is  basically  a  set  of  simultaneous  quadratic  algebraic 
equations.  If  the  variables  are  regarded  as  being  the  relative  populations 
x.  *  n./n,  where  n  is  the  total  CO  population,  then  for  an  M-level  approxi¬ 
mation  there  will  be  a  system  of  (M+l)  equations , 


Vx>  =  ?A2'jkXjXk+  ?B2,kXk'° 

*«»>,  EA.  X.X  +2£b.  „ 

i  jk  i,  jk  j  k  k  i,  k  k 


,  x  =  0 
k  k 


(X  )  =  ^Aw  ..  x.x  +  S  Bw  .  x  =  0 
L  jk  M,jk  ]  k  k  M,  k  k 


(3.75) 


fM+l  (  x  )  =  +  Xj  +  x2  +  *  *  *  _1  =  ° 
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where  f  ,  f  ,  ...  f  are  the  rates  for  levels  1,  2,  ....  M,  and  the  (M+  1)  st 

1  2  M  rp 

equation  expresses  total  particle  conservation,  x  =  (x  ,  x  ,  ....  x  ,  x  ) 

|v  “  ^  ®  IJI 

is  an  (M+l)  dimension  column  vector,  and  if  a  vector  7  =  (f  ,  f^, . . .  f^,  fj^+j) 
is  defined,  the  molecular  kinetic  equation  can  be  written  concisely  as 


f  (5)  =  0. 


(3.75a) 


The  technique  that  was  employed  for  solving  this  system  of  equations 

68 

numerically  is  a  generalisation  of  the  Newton -Raph son  method  for  solution 

of  nonlinear  equations.  Define  the  (M+l)  x  (M+l)  Jacobian  matrix 


1^.(3)=  dfjf*)/ dx.. 


(3.76) 


If  x  is  a  guess  for  the  solution  vector  a  which  satisfies  f  (a)  =  0,  then 
a  "better"  guess  is  given  by  the  vector  x*,  defined  by 


■*  1  -o  _/*o -1  . -o. 

x  =  x  -  J(x  )  *  f(x  ) 


(3.77) 


Starting  with  some  initial  guess  vector  x  ,  a  sequence  of  vectors 

x  ,  i  ,  ,,,  x  is  generated  in  this  way,  and  successive  iterations  will 

usually  converage  to  the  solution  "$  ,  generally  with  very  minimal  restrictive 

conditions.  For  the  present  set  of  equations  with  M  =  40,  this  procedure 

typically  required  ~10  iterations  to  produce  a  solution  accurate  to  one  part 
7 

in  10  for  all  of  the  populations,  and  computation  time  on  a  CDC  6600  was 
about  60  sec.  These  characteristics  correspond  to  very  simple  forms  for 
the  initial  guess  ^x0  (e.  g. ,  "x°  =  (0,  0,  0, ...  0, 1)^),  and  could  probably 
be  improved  with  a  better  starting  point. 
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3.6.4  Transient  Analysis.  For  situations  where  the  steady  state  has  not 
been  attained,  it  is  necessary  to  solve  the  simultaneous  set  of  nonlinear 
first  order  differential  equations,  Eq.  (3.1a),  subject  to  a  set  of  prescribed 
initial  conditions.  In  that  case,  the  theoretical  analysis  described  in  the 
preceding  sections  would  have  to  be  extended.  However,  the  mods  of 
operation  for  which  such  transient  effects  would  be  expected  to  play  an 
important  role  would  be  a  short  pulse  regime,  where  the  time  scale  is 
typically  less  than  (or  comparable  to)  the  fastest  relaxation  lifetimes  that 
dominate  the  population  distributions.  For  high  pressure  molecular 
gases,  the  fastest  mechanisms  that  determine  that  distribution  are  the 
near -resonant  YVT  exchange  collisions,  with  cross -relaxation  rates 
typically  of  the  order  of  pr  —400  0  sector  r.  For  a  CO  partial  pressure 
p__ ,  it  would  be  expected  that  a  transient  analysis  based  on  the  set  of 

VV 

differential  equations  would  only  be  required  for  pulsed  operation  with  times 
scales  £  (pr)/p£Q.  Therefore,  for  long  pulse  or  cw  operation,  it  is 
to  be  expected  that  temporal  effects  could  be  treated  in  a  quasi-static 
approximation,  where  parameters  which  are  slowly  varying  functions  of 
time  (e.g. ,  the  kinetic  temperature  T(t)  varies  during  a  pulse)  can  be 
regarded  as  fixed  values  which  determine,  at  any  given  time  t,  an 
instantaneously  established  steady-state  population  distribution.  Thus, 
in  the  same  spirit,  quantities  such  as  gain  and  saturation  intensities 
will  be  quasi-steady,  slowly  varying  functions  of  time  during  the  pulse. 
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4.0  PLASMA  CHARACTERISTICS  AND  ELECTRON  KINETICS 

The  kinetic  model  described  for  calculation  of  the  steady-state  vibrational 
populations  has  assumed,  for  the  present,  that  the  distribution  of  electrons 
in  the  medium  is  a  Boltzmann  function  at  temperature  Te.  In  actual  practice, 
this  distribution  can  depart  significantly  from  the  Maxwell-Boltzmann  form, 
and  it  would  be  an  improvement  to  calculate  its  energy  dependence  more 
accurately.  For  the  weakly  ionized  plasmas  that  are  characteristic  of 
electrical  discharge  laser  systems  (i.  e. ,  electron  densities  much  less 
that  neutral  molecular  densities)  it  is  the  inelastic  vibrational  and  electronic 
excitation  processes  that  dominate  the  kinetics  of  the  electron  gas.  The 
dominant  mechanism  that  determines  the  form  of  the  electron  distribution 
is  collisions  between  electrons  and  molecules  that  result  in  vibrational  or 
electronic  state  excitation.  The  effect  of  these  inelastic  collisions  is  to 
extract  energy  from  the  electrons  as  they  are  accelerated  across  the  region 
of  the  applied  electric  field.  The  combined  effect  of  the  electric  field  and  the 
inelastic  collisions  is  to  limit  the  electron  drift  velocity,  which  is  determined 
by  the  value  of  (E/N)  and  the  molecular  gas  composition.  Since  only  low- 
energy  electrons  (~  1-3  eV)  are  effective  for  vibrational  excitation,  it  is 
desirable  to  tailor  the  (E/N)  and  the  gas  composition  to  optimally  achieve 
this  objective. 

The  structure  of  the  electron  distribution  will  typically  reflect  the  structure 
of  the  electron -molecule  energy  exchange  processes,  and  thus,  it  will 
depart  most  significantly  from  a  Boltzmann  in  the  energy  regions  where 
inelastic  scattering  cross-sections  are  appreciable.  For  low  density 
plasmas  with  mean  electron  energies  of  a  few  eV  or  less,  elastic  electron 
scattering  and  inelastic  electron-molecule  rotational  excitation  have  a 
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negligible  effect  on  the  determination  of  the  distribution.  Since  the 
electrons  for  such  low  energy,  weakly  ionized  plasmas  have  no  mechanism 
for  complete  thermalization,  a  departure  from  a  Boltzmann  distribution 
is  to  be  expected. 

For  the  present  situation,  the  electron  kinetic  problem  involves  the  solution 
of  the  Boltzmann  equation^’  ^for  the  energy  distribution  function  that 
results  for  a  given  gas  composition  in  the  presence  of  an  applied  electric 
field  E.  In  general,  the  solution  of  the  Boltzmann  equation  is  often  a  practical 
impossibility,  since  complete  information  about  all  of  the  collision  processes 
is  usually  unavailable.  However,  for  the  gas  constituents  (e.g.  CO,  ^,002* 
He,  Ar,  Xe)  and  electron  energy  ranges  (~  few  eV)  that  are  typical  for  mole¬ 
cular  laser  systems,  reliable  cross-sections  for  the  important  scattering 
processes  have  been  measured  experimentally.^  ^  Nighan  et  al  ^ 
have  carried  out  calculations  for  the  electron  energy  distribution  functions 
for  various  gas  mixtures,  based  on  the  assumptions  that  the  E/N  of  the  plasma 
is  low  enough  that  the  average  electron  energy  is  only  a  few  eV,  and  that 
the  electron  density  is  low  enough  that  the  approximation  of  a  weakly 
ionized  plasma  is  valid.  Figure  4. 1  to  4.3  ,  taken  from  Reference  71,  74 
illustrate  the  significant  departure  from  Boltzmann  form  that  can  result 
for  the  electron  energy  distribution  under  various  conditions. 

Knowledge  of  the  electron  energy  distribution  as  a  function  of  E/N  and  gas 
composition  makes  it  possible  to  calculate  the  fractional  amounts  of  electrical 
energy  input  that  are  deposited  into  each  of  the  possible  excitation  modes 
present.  The  fractional  amount  of  energy  transferred  to  molecular  vibrational 
states,  in  conjunction  with  the  quantum  extraction  efficiency  discussed  pre¬ 
viously,  will  then  determine  the  total  efficiency  of  conversion  of  electrical 
input  to  optical  output  power.  Clearly,  it  is  desirable  to  operate  with 
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Figure  4. 1.  Calculated  electron  energy  distributions  for  a 
CC^-^-He  gas  mixture  (1-1-8),  taken  from 
Reference  71. 
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Figure  4.2.  Comparison  of  calculated  electron  energy  distribution 

functions  with  a  Maxwellian  for  Ng,  CO,  and  CO2  (taken 
from  Reference  71).  This  exhibits  the  significant  departure 
from  a  Boltzmann  distribution  that  can  occur  for  a  weakly 
ionized  plasma. 
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Figure  4.  3.  Electron  energy  distribution  functions  for  CO,  CO 

7A  ^ various  E/N  values,  taken  from  Reference 

functiN  Ul/2f(u)  corresponds  to  the  distribution 

function  as  defined  in  this  report). 
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discharge  characteristics  that  will  optimize  the  transfer  of  electrical 
energy  to  the  vibrational  levels,  rather  than  to  electronic  or  ionized 
states.  This  is  a  difficult  objective  to  achieve  for  the  conventional 
longitudinal  electric  discharge  configurations  that  are  usually  employed 
for  conventional  CO  and  CO2  lasers,  since  the  electron  distribution 
function  must  adjust  itself  in  such  a  way  that  average  electron  energies 
are  high  enough  to  provide  the  ionization  that  is  required  to  sustain  the 
discharge.  Thus,  one  of  the  ways  that  is  often  successfully  employed 
for  enhancing  gas  laser  efficiencies  is  to  optimize  the  mixture  with  gases 
that  have  a  low  ionization  potential,  such  as  Xe  or  Hg  vapor  ’  . 

Although  plasma  control  by  optimization  of  gas  constitutents  is  a  useful 
technique  for  improving  laser  operating  characteristics,  the  degree  of 
control  that  this  can  provide  is  usually  very  limited. 

The  optimum  gas  mixture  for  conventional  CO  lasers  contains  several  other 

constituents  in  addition  to  CO.  In  a  slowly  flowing  system  He  is  well  known 

to  be  beneficial  for  cooling  the  molecular  kinetic  temperature.  The  effect  of 

Xe  is  to  reduce  the  average  electron  energy  by  lowering  the  high  energy 

tail  of  the  electron  energy  distribution.  Use  of  Xe  is,  therefore,  a  crude 

method  for  altering  the  E/N  of  the  discharge.  Since  Xe  has  the  lowest 

ionization  potential  of  any  of  the  constituents  typically  used  in  a  CO  laser 

mixture,  it  may  also  be  useful  for  preventing  CO  decomposition  due  to 

collisions  with  high  energy  electrons.  According  to  the  measured 

29 

cross -sections  of  Shulz  (Figure  3.1  ),  electron  excitation  of  CO  is 
effective  only  for  £  8  vibrational  quanta.  Qualitatively,  the  cross-sections 
for  excitation  of  Ng  are  quite  similar  to  those  for  CO,  although  they  are 
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smaller  by  about  a  factor  of  six.  Hence,  would  be  expected  to  have  an 
effect  on  the  electron  distribution  qualitatively  similar  (in  energy  dependence) 
to  that  for  CO.  Therefore,  addition  of  Ng  to  a  CO  laser  mixture  can  be  of 
benefit,  since  for  the  proper  values  of  E/N,  it  would  be  as  efficient  for  stor¬ 
ing  electrical  energy  in  vibrational  states  as  is  CO.  Since  the  levels  of  N£  are 
nearly  resonant  with  those  for  CO,  this  vibrational  energy  can  be  transferred 
by  V-V-T  cross -relaxation  to  the  CO  levels  and  would  not  be  last  (since  V-T 
relaxation  is  much  slower).  Two  important  factors  to  realize,  however,  are 
that  N2  can  be  added  to  alter  the  E/N  of  the  plasma  without  loss  of  electron- 
vibration  energy  transfer  efficiency,  and  its  addition  does  not  directly  affect 
the  CO-CO  V-V-T  anharmonic  pumping  rates,  which  would  have  occured  if  we 
had  merely  increased  the  CO  partial  pressure.  Since  V-V-T  pumping  in  the  N2- 
N2  system  is  much  slower  than  that  for  N2-CO,  most  of  the  transfer  from  N2 
will  be  from  the  lowest  level.  The  beneficial  effects  that  result  from  the 
simultaneous  pressure  of  N2  and  Xe  may  be  due  to  the  fact  that  and  Xe 
alter  the  plasma  characteristics  in  qualitatively  different  ways.  For  example, 
the  fractional  conversion  of  electrical  to  vibrational  energy  is  improved  not 
only  upon  lowering  the  E/N,  but  is  evidently  enhanced  also  by  an  optimum 
gas  composition,  since  the  performance  can  be  improved  by  the  simultaneous 
presence  of  both  of  these  constituents.  An  equivalent  enhancement  cannot  be 
obtained  merely  by  increasing  the  fractional  amount  of  Xe  in  the  gas  mixture. 
Since  the  addition  of  N2  does  not  cause  degradation  of  the  fractional  transfer 
of  energy  to  vibrations,  it  may  have  advantages  that  an  addition  of  Xe  cannot 
provide.  For  the  plasma  characteristics,  the  effect  of  the  addition  of  N2 
could  probably  be  accomplished  equivalently  by  increasing  the  partial 
pressure  of  CO  (although  the  average  electron  energy  would  not  be  the  same 
function  of  E/N).  However,  the  optimum  concentration  of  CO  in  the  laser 
mixture  is  also  determined  by  other  constraints,  such  as  the  rate  of  V-V-T 
pumping,  V-T  relaxation  decomposition,  etc.  Thus,  the  optimum  mixture 
seems  to  require  several  ingredients  (typically  CO,  N2,  He  and  Xe).  The 
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addition  of  O ^  often  plays  a  very  sensitive  role,  inasmuch  as  very  small 
amounts  can  have  a  critical  effect  on  the  operation.  A  small  amount  of  O2 
is  helpful,  because  it  probably  minimizes  CO  decomposition  by  driving  the 
reaction  COi?"  C  +  O  to  the  left  and  prohibits  carbon  deposits.  Beyond  a 
critical  amount,  the  presence  of  O2  can  have  a  deleterious  effect  on  the  laser 
performance.  This  may  be  due  to  a  drastic  alteration  of  the  plasma  character¬ 
istics  (since  O2  has  a  large  electron  attachment  coefficient  of  magnitude 
greater  than  that  for  most  gases),  or  it  may  result  in  the  formation  of  CO2 
(at  least  for  room  temperature  operation).  Since  such  minute  amounts  can 

have  such  an  important  effect,  the  former  explanation  is  probably  correct. 

71-74 

Nighan  has  shown  that  for  electron  energies  lower  than  those  necessary 

for  a  self-sustaining  glow  discharge,  the  efficiency  of  electrical  energy 
transfer  to  the  vibrational  levels  can  be  greatly  enhanced,  and  his  analyses 
show  the  extreme  sensitivity  in  some  cases  of  the  vibrational  excitation  rates 
upon  plasma  characteristics.  Figure  4.4  and  4.  5,  taken  from  Reference  71“ 
74  illustrate  these  results.  Clearly,  in  order  to  best  achieve  an  optimiza¬ 
tion  of  the  plasma  characteristics  of  a  laser  system,  it  would  be  desirable 
to  have  an  independent  control  over  the  ionization  process  and  the  electrical 
excitation  of  the  molecular  vibrational  levels,  since  then  the  E/N  could  be 
adjusted  to  have  the  most  favorable  value.  One  sophisticated  technique  for  doing 
that  is  to  obtain  the  initial  pre-ionization  with  a  high  energy  electron  beam, 
which  creates  a  volume  ionized  gas  of  low  energy  secondary  electrons  which 
are  then  accelerated  by  an  applied  electric  field  to  sustain  the  discharge 
current.  For  the  case  of  CO,  Figure  4.  5  indicates  that  optimization 
of  energy  transfer  to  vibrational  levels  is  not  a  sensitive  function  of  E/N 
below  a  certain  minimum  value. 

For  the  case  of  CO2  laser  systems,  the  molecular  kinetics  are  much 
simpler  than  for  the  case  of  CO,  where  many  vibrational  levels  must  be 
considered.  The  electron  distribution  function  must  depend  upon  the 
populations  of  excited  vibrational  states  of  the  molecular  species  if  there 
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Figure  4.4.  Fractional  electron  power  transfer  to  competing 
excitation  processes  for  electrically  excited  CO2- 
N2~He  laser  mixture,  taken  from  Reference  71. 
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5.  Fractional  electron  power  transferred  to 
vioratxonal  and  electronic  state  excitation 
for  (a)  N,,  and  (b)  CO  (taken  from  Reference 
74). 


exist  inelastic  excitation  processes  from  such  states  that  are  important  for 
•  extracting  energy  from  the  electron  gas.  Thus,  a  complete  model  for  an 

electrical  discharge  CO  laser  would  require  a  coupled  set  of  equations  for 
the  molecular  and  electron  kinetics.  Our  experiments  have  shown  that  the 
presence  or  absence  of  lasing  may  significantly  change  the  electron  distribu¬ 
tion  function;  this  was  demonstrated  by  the  fact  that  the  visible  CO  emission 
from  upper  electronic  states  sharply  increases  when  the  laser  cavity  is 
blocked,  suggesting  that  the  average  electron  energy  has  been  increased  by 
the  altered  population  distribution  of  the  CO  vibrational  states.  Since  the 
vibrational  population  distribution  is  different  when  lasing  is  prohibited,  the 

relative  importance  of  the  vibrational  versus  electron  excitation  mechanisms 

74 

for  the  electrons  would  be  changed.  Nighan  has,  in  fact,  carried  out  cal¬ 
culations  for  the  electron  distribution  that  arises  in  vibrationally  excited  N^, 
assuming  that  electron  excitation  collisions  with  N^(v  =  1)  are  also  included. 
Fig.  4.6,  taken  from  Ref.  74,  illustrates  the  effect  that  the  vibrational 
population  distribution  can  have  on  the  electron  kinetics.  Thus,  for  the  CO 
laser,  it  is  probably  not  completely  accurate  to  characterize  the  plasma  only 
by  the  gas  composition  and  the  value  of  E/N;  the  details  of  the  vibrational 
population  distribution  may  also  be  important  for  the  coupled  system  of  mole¬ 
cules  and  electrons.  For  the  case  of  the  CO^  laser  system,  extremely  rapid 
cross*- relaxation  processes  make  it  possible  to  describe  the  molecular  kine¬ 
tics  in  terms  of  only  a  very  few  molecular  vibrational  levels,  so  that  the  sig¬ 
nificant  excitation  processes  that  determine  the  electron  distribution  occur 
mainly  from  the  ground  state. 

In  addition  to  further  improvements  on  the  molecular  kinetic  model  for  the 
CO  system,  it  would  be  desirable  to  undertake  theoretical  calculations  of  the 
electron  energy  distribution  functions.  This  would  be  especially  valuable  for 
treating  the  combined  system  of  molecular  and  electron  kinetics  on  a  con¬ 
sistent  and  equal  basis,  and  is  necessary  for  constructing  a  reliable  model 
for  electric  discharge  CO  laser  systems. 
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.6  .  Electron  energy  distribution  functions  in 
Ng  calculated  for  various  vibrational 
temperatures  (taken  from  Reference  74). 
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5 .  0  SUMMARY  AND  C ONC LUSIONS 


A  complete  theoretical  understanding  of  the  operating  characteristics  of  an 
electrically  excited  CO  laser  system  requires  a  detailed  model  for  the  kine¬ 
tics  of  both  the  electrons  and  the  molecular  vibrational  levels.  Although  a 
completely  .self-consistent  analysis  should  attempt  to  treat  these  two  coupled 
systems  on  an  equal  basis,  there  is  much  useful  information  that  can  be  ob¬ 
tained  by  analyzing  them  independently.  The  basic  reason  for  this  is  that 
the  characteristics  of  one  system  are  only  reflected  somewhat  insensitively 
in  the  resulting  properties  for  the  other.  For  example,  the  detailed  struc¬ 
ture  of  the  electron  energy  distribution  does  not  critically  affect  the  integra¬ 
tions  over  the  cross  sections  that  give  the  rates  for  electron  pumping  of  the 
molecular  vibrational  levels.  That  is,  even  though  the  electron  energy  dis¬ 
tribution  function  may  display  very  important  behaviour  in  the  energy  regions 
of  interest  for  vibrational  excitation,  use  of  an  (incorrect)  Boltzmann  distri¬ 
bution  can  give  approximately  the  correct  rates  if  the  electron  temperature 
is  selected  properly.  On  the  other  hand,  the  detailed  structure  of  the  vibra¬ 
tional  population  distribution  does  not  have  a  profound  effect  upon  the  energy 
dependence  of  the  electron  distribution,  although  experimental  evidence  indi¬ 
cates  that  the  plasma  characteristics  may  change  in  the  presence  or  absence 
of  lasing.  The  tacit  assumption  for  an  independent  calculation  for  the  elec¬ 
tron  distribution  is  that  most  of  the  molecules  are  in  the  ground  state,  and 
that  the  significant  processes  that  determine  the  energy  dependence  are  in¬ 
elastic  excitations  from  the  ground  state.  To  the  extent  that  all  of  these 
approximations  are  valid,  a  reasonable  theoretical  picture  of  electrical  CO 
laser  systems  can  be  obtained  from  results  obtained  using  independent  kine¬ 
tic  models  for  the  electrons  and  the  molecular  vibrational  levels. 

As  far  as  the  molecular  kinetics  are  concerned,  there  are  several  processes 
of  pumping,  relaxation,  and  cross -relaxation  that  must  be  included  to  ex¬ 
plain  theoretically  the  population  distributions  that  are  produced.  Although 
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the  initial  source  of  energy  input  is  the  electrical  excitation  of  the  lowest  vi¬ 
brational  levels,  the  dominant  vibration-vibration  exchange  collisions  (i.  e.  , 
the  "anharmonic  pumping"  mechanism)  are  responsible  for  rapid  redistribu¬ 
tion  of  this  energy  to  higher  levels.  In  this  report,  a  discussion  of  several 
important  aspects  of  the  CO  laser  has  been  given,  along  with  a  description 
of  the  molecular  kinotic  model  and  computer  program  that  has  been  constructed 
at  NCL  for  calculations  of  the  characteristics  of  a  steady-state  CO  laser  am¬ 
plifier  system.  Although  there  presently  remains  a  need  for  refining  some 
of  the  experimental  constants  and  theories  of  rates  needed  for  this  model, 
the  program  is  essentially  complete.  From  a  reasonably  general  input  de¬ 
scription  of  the  experimental  conditions,  it  is  designed  to  calculate  all  of  the 
information  about  the  vibrational  population  distribution  that  is  necessary  to 
completely  describe  the  radiative  characteristics.  Small-signal  gain  meas¬ 
urements  as  well  as  measurements  of  relaxation  rates  using  double  Q-switch- 

21 

ing  and  double  probe  techniques  are  planned  for  the  near  future  ,  and  should 
help  to  finalize  some  of  the  constants  required  for  construction  of  a  reliable 
model. 

For  completeness,  a  discussion  of  several  important  aspects  of  the  CO  laser 
has  been  given.  In  addition  to  a  qualitative  and  quantitative  discussion  of  the 
molecular  and  electron  kinetics,  with  formulation  and  numerical  solution  of 
the  master  equation  for  the  steady  state  vibrational  populations,  several  use¬ 
ful  equations  for  the  spectroscopic  and  radiative  properties  have  also  been 
included.  The  purpose  for  expanding  the  scope  of  this  report  somewhat  be¬ 
yond  the  specific  topic  of  molecular  kinetics  has  been  to  present  a  unified 
treatment  of  the  physics  of  the  electrically  excited  CO  laser,  as  well  as  to 
provide  a  convenient  summary  of  several  relevant  equations  and  references. 
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